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ABSTRACT 


The last week of April and first week of May 1988 along the central California coast 
were characterized as a period of strong coastal upwelling produced by moderate to 
strong northwesterly winds present throughout the period. A product of this upwelling 
event was the manifestation of a southward geostrophic jet which extended to a distance 
of approximately 50 km from the coast. 

During the period from 08 to 11 May 1988, a hydrographic survey consisting of 17 
stations each occupied four to five times was conducted within the Monterey Bay. 
Internal waves, with amplitudes of up to 30 m were present throughout the period and 
effectively masked the mean signal, implying that averaging is essential to avoid aliasing. 
The CTD data were averaged to estimate the mean field during this time frame. ADCP 
data, also acquired during this period, were also averaged. The mean (low field and dy¬ 
namic topography implied anticyclonic surface flow with cyclonic flow at 200 m depth, 
both roughly centered over the axis of the Monterey Submarine Canyon near the mouth 
of the Bay. Flow at depth appeared to be Canyon "trapped." Surface current speeds 
were on the order of 10 to 15 cm s -1 with somewhat slower speeds at depth of approxi¬ 
mately 10 cm srK ADCP derived mean flows compared favorably with geostrophic mean 
flows in all areas except one, the deep outflow region along the northern wall of the 
Canyon. Calculated standard errors were smaller than the mean signal for both data 
sets. The smaller standard errors calculated for the ADCP data implied that this method 
resolved the mean signal better than did the geostrophic calculations. 

Application of the geostrophic mean field to the open boundary of a two layer, 
primitive equation numerical ocean model yielded flows similar to those described above. 
The inclusion of linear bottom friction was shown to be particularly important in limit¬ 
ing deep flow to the confines of the Canyon. Wind stress forcing experiments indicated 
that a strong wind field may influence surface circulation in the Bay. 

Mean flow fields compared favorably with results of Klinck's (1989) three-level 
model of geostrophic adjustment of flow over a "narrow" submarine canyon. Inter¬ 
actions between the coastal upwelling geostrophic jet and the Monterey Submarine 
Canyon is believed to have been a major mechanism responsible for producing the ob¬ 
served mean flow. 
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I. INTRODUCTION 


A. THE MONTEREY BAY AND VICINITY 

1. Dimensions 

The Monterey Bay, located on the central California coast, is basically orien¬ 
tated north-south with the city of Santa Cruz located on the northern shore, Moss 
Landing to the east, and the Monterey Peninsula situated on the southern shore. To the 
west, the bay is open to the sea. A line drawn between Santa Cruz and the Monterey 
Peninsula measures 37 km, and Moss Landing lies approximately 20 km east of this line. 
The bay contains two bights (Breaker and Broenkow, 1989), one situated along the 
northern coast between Santa Cruz and Aptos, the other between the Monterey Penin¬ 
sula and Fort Ord (Figure 1 on page 2). The area of the bay is approximately 550 km 2 
constituting a volume of approximately 420 km 3 (Martin, 1964). 

Monterey Bay is essentially a shallow water bay (less than 100 m deep) split into 
northern and southern shelf regions by the Monterey Submarine Canyon. This canyon 
is the largest bathymetric feature on the west coast, comparable in shape and relief to 
the Grand Canyon of the Colorado (Shepard, 1966). The canyon's axis runs roughly 
east-west within the bay, originating near the surface at the canyon head near Moss 
Landing where it is 5 km in width, then widens to 15 km, rapidly deepens to 1500 to 2000 
m and begins to meander past the mouth of the bay, northwest of the Monterey Penin¬ 
sula. The two largest tributaries entering the Monterey Canyon are the Soquel and 
Carmel Canyons, the former intersecting the axis near the center of the bay from the 
northeast and the latter entering it from the south just seaward of the Monterey Penin¬ 
sula. Past the entrance to the bay, the axis runs southwest for 20 km, then heads 
seaward around a large fan feature to the north and continues deepening to over 3000 
m. The relief of this feature is truly dramatic (Figure 2 on page 3). 

2. Currents of the Central California Coast 

Because circulation within the bay »s undoubtably affected by the flow offshore 
(Breaker and Broenkow, 1989), a brief description of major current systems in the vi¬ 
cinity of the California coast, and their generation mechanisms, is provided. The three 
major currents influencing the California coast arc the California Current (CC), the 
California Undercurrent (CUC), and the Davidson Current. These names are usually 
applied to features existing in the long term or climatological sense (Hickey, 1979). 
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Figure 1. Hie Monterey Bay (Source: Breaker and Broenkow, 1989). 
a. The California Current 

'I he California Cuirent refers to the broad, diffuse, cquatorward wind driven 
flow of the eastern boundary of the North Pacific Subtropical anticylonic gyre. This 
current dominates (in the mean) the eastern Pacific How and extends south from the 
North Pacific to approximately 25’ N, where it turns to the southwest. This relatively 
fresh (Pacific Subarctic Watci) cuncnt is usually strongest within a few hundred kilo- 
mctcis of the coast. The western boundary, of inteicst in this study, has arbitrarily been 
set at 1000 km from shore. Mean velocities range fiom 1 to 10 cm j 1 , typical for most 
eastern boundary currents. Of more interest to this study aic the CUC and the 
Davidson current, both located in the near shore legions, typically cast of the CC 
(Hickey, 1979). 










Figure 2. Bottom topography of the Monterey Bay Submarine Canyon and 
vicinity: lop figure extends seaward to 112 km, bottom to 60 km (fig¬ 


ures derived from NOAA/NOS 250 m resolution gridded data). 
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b. The California Undercurrent and Davidson Current 

The CUC is the mean poleward flowing current with mean core velocities 
of 10 cm s -1 that is typically found along the continental slope at depths ranging from 
100 to 300 m, The Davidson Current refers to the northward flowing surface current 
found along the coast with velocities on the order of 10 cm s~ l and a width of approxi¬ 
mately 80 km (McCreary et al., 1987). 

Prevailing winds out of the northwest, generated by a quasi-stationary at¬ 
mospheric high pressure system (Figure 3 on page 5) over the north central Pacific, are 
present year round over the California coast. The North Pacific High typically 
strengthens during the summer months while the thermal low pressure over the 
California coast, particularly in the central valley, deepens as the land warms. The 
combined effects of these phenomena results in a tighter atmospheric pressure gradient 
producing strong northwesterly winds along the California coast during the summer 
months (Huyer, 1983). Halliwell and Allen (1987) observed that high pressure frequently 
builds northeast across California after the passage of an atmospheric cyclone and its 
associated cold front generating periods of strong northwest winds and upwelling events 
along the coast. They also noted that wind maxima tended to propagate northward 
along the coast. Chelton (1984) noted that the spatial structure of the wind stress field 
is highly coherent in the alongshore direction and that the magnitude varies seasonally 
with the maximum intensity in May (Figure 4 on page 7). Since the strongest compo¬ 
nents are found approximately 200 km offshore, a positive wind stress curl (curl r) is 
present year round in the nearshore regions (Nelson, 1977). These winds, through the 
Ekman transport, produce near shore upwelling resulting in dynamic height fields which 
slope downward towards the coast. The result is a geostrophically balanced surface jet 
which flows in the direction of the wind (i.e. equatorward) and a poleward undercurrent 
(the CUC) at approximately 100 to 300 m in depth. Hickey (1979) observed that this 
surface flow is strongest along the Pacific northwest coast during the spring with maxi¬ 
mum flow typically located within the internal Rossby radius of deformation, between 
5 and 25 km from the coast, and noted that variations in width, location, and strength 
can occur quite rapidly (on the order of days) in response to local changes in wind stress 
(t). She also observed that during periods of weak northerly or southerly winds, this 
surface jet either weakens or disappears altogether. Hickey (1979) finally noted that the 
dynamics of the nearshore southward flow observed ofT the California coast are likely 
similar to the flows noted along the Pacific northwest coast, 
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Figure 3. Long-term mean sea level atmospheric pressure over the no.th 
Pacific: Top figure is for January, bottom figure is for July (Courcc: 
Iluyci, 1983). 


During the summer, observed current 1 : along the coast arc in general 
agreement with i'kman transpoit theory but during the winter poleward flow (the 
Davidson Current) occurs along the coast, against the surface wind sticss (McCreary ct 
al., 1987). McCreaiy et al. (1987) cite four possible mechanisms, along with pros and 
cons of each, which might explain the appearance of the Davidson curtcnt, namely; 

• Poleward propagation of coastally trapped waves. 

• Relaxation of local winds. Sustained northwesterly winds produce a poleward 
pressure gradient. When these winds relax, a poleward cuircnt will be driven by 
this pressure gradient. 'I his theory is consistent t’ c appearance of the Davidson 
Current during the winter, when mean equatorwaid suifacc wind stress decreases. 
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• Positive wind stress curl. Sverdrup (1947) surmised that the vertically integrated 
transport in the meridional direction (M y ) is proportional to the vertical component 
of the curl of the wind stress, or: 
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where fi = djldy is the meridional change in the Coriolis parameter (/), ^ is the 
stream function (d^/dx > 0 implies northward flow), and t x and r y are the x and y 
components of the wind stress respectively. Thus, positive curl t should drive a 
poleward current, according to Sverdrup theory. The problem with this mechanism 
however, is the fact the Davidson Current is usually found during the winter 
months, a time period when curl t is typically the weakest. 

• Thermohaline forcing. The north-south density gradient between subtropical and 
polar regions creates lower dynamic heights over the poles and relatively higher 
dynamic heights in the equatorial region. This situation drives a mid-ocean east¬ 
ward geostrophic current which should turn cyclonically (poleward) when it en¬ 
counters an eastern boundary such as the California coast. They argue that this 
mechanism is likely to be responsible for driving the Leeuwin Current off western 
Australia; however, since wintertime sea-surface topography off tilt California 
coast differs significantly from that observed off the west Australian coast, (i.e. 
dynamic heights bend toward the equator vice the pole off California), this mech¬ 
anism appears to play a role which is subordinate to wind forcing in this area. 


Using both wind stress curl and curl free wind stress as inputs into a coastal model, 
McC r eary, et al. (1987), concluded that positive wind stress curl was the most likely 
mechanism responsible for driving the Davidson Current during the winter. During the 
summer months the along shore wind stress (r y ) is strong enough to negate the effects 
of Sverdrup transport and therefore dominates the surface flow' field. Hickey (1979) 
states that oceanic responses to sign changes in curl t w'ould be slow, i.e. on the order 
of two to three months. 

Chelton (1984), through the averaging of 23 years of hydrographic data ac¬ 
quired through the CalCOFI program, found the surface flow field to be closely coupled 
to the w'ind field. Analysis of the derived dynamic height fields at the surface and at 150 
dbar relative to 500 dbar at Monterey (station 67.57) showed surface flow to be south¬ 
ward from February to June and from August to November with polew'ard flow from 
November to February and again in June to July. Subsurface flow (CUC) towards the 
north w’as present year round except for weak southward flow in April and May. Further 
south, at Point Conception, no reversal in the CUC w'as noted, and thus Chelton con¬ 
cluded that subsurface convergence must be taking place somewhere between Pt Con- 
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Figure 4. Long term average wind stress (top) and wind stress curl (bottom) along 
the California coast: Averages were obtained from 100 years of ship 
observations. The shaded region corresponds to positive wind stress curl 
(Source: Nelson, 1977). 

ccption and Pt. Sur, although the CalCOFI grid spacing was too coarse to resolve this 
convergence. In general, Chelton (1984) concludes that his findings support the general 
description of an equatorward CC, a poleward flowing CUC, and a coastal surface cur¬ 
rent whose direction of flow is seasonally dependent. 

In summary, since coastal circulation undoubtably influences, and possibly 
drives, flow within the Monterey Bay, full knowledge of the near shore processes which 
take place along the central California coast is crucial prior to attempting to understand 
(low patterns within the bay itself. For instance, surface flow within the Monterey Bay 
during winter, when the Davidson Current is present, more than likely differs from flow 
patterns within the bay during spring upweiling periods when the southward geostrophic 
jet is present. Also, deep flow within the Monterey Bay Submarine Canyon is probably 
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influenced to some extent by the CUC, although the effects of this current on circulation 
within the bay does not appear to be very well understood at present. 

B. SYNOPSIS OF PREVIOUS MONTEREY BAY STUDIES 
1. Direct Observations 

Numerous studies have been conducted to ascertain the characteristics of the 
circulation within Monterey Bay utilizing a wide variety of methods including current 
meters, drift bottles, drogues, hydrographic data, thermistors, satellite imagery, aircraft 
observations, and sediment and nutrient distributions. Breaker and Broenkow (1989) 
provide an excellent in-depth summary of past research conducted in the bay, including 
comments relating these works to their own findings. Highlights of some of these major 
studies are provided. 

a. Early Studies 

One of the earliest studies of flow patterns within Monterey Bay was con¬ 
ducted by Bigelow and Leslie (1930) using hydrographic data obtained throughout the 
bay over the period 30 June to 24 July 1928. They observed higher sea surface temper¬ 
atures over the shelf areas of the bay and somewhat depressed temperatures over the 
canyon implying that local heating played a major role during the summer months. 
Analysis of the cross sectional temperature field also led them to believe that upwelling 
from within the canyon occurred and subsequently spread out into the shelf regions of 
the bay. Dynamic height analysis (0 dbar referenced to 500 dbar) revealed the presence 
of high dynamic heights near the center of the canyon implying anticyclonic flow. They 
concluded that circulation within Monterey Bay varies greatly and is dramatically influ¬ 
enced by week to week changes in small scale local events (heating, land-sea breeze, 
winds, canyon upwelling, etc.). 

Skogsberg (1936), through analysis of hydrographic data acquired between 
1929 and 1933 in the southern portion of Monterey Bay, concluded that three separate 
cycles occurred annually. The first he described as a "cold water" or upwelling phase 
which lasted from February through November, a "warm water phase" attributed to light 
winds and onshore movement of warm oceanic waters from August through November, 
and finally, a "low thermal gradient phase" which lasted from December to February and 
coincided with the manifestation of the northward flowing Davidson Current. He em¬ 
phasized the highly variable nature of flow vithin the bay stating instances where coastal 
flow was northward through the bay and other instances where it was southward. 
Through these thermal studies, Skogsberg finally noted that temperature trends at lower 
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levels were quite often opposite to those encountered at shallower depths suggesting that 
flows may, at times, be independent and opposite at separate levels within the bay. 
b. Surface Current Studies 

Utilizing nutrient displacement and mass fields, Lasley (1977) concluded 
that surface flow within the bay is predominantly to the north, frequently entering the 
bay around Pt. Pinos and exiting near Pt. Santa Cruz. Utilizing property distributions 
assimilated over a 27-month time span, Smethie (1973), and Broenkow and Smethie 
(1978), also inferred surface flow within the bay to be northward the majority of the 
time, in agreement with Lasley. They expected northward flow during the fall and winter; 
however, in light of the northwesterly winds present during the spring and summer, 
prevailing northward surface flow within the Monterey Bay was contrary to their ex¬ 
pectations during this time frame. Current meter data acquired by Engineering Science, 
Inc. (ESI, 1978) in the area southwest of the Salinas River at depths of 9 and 30 m over 
the period from January 1976 to January 1977, indicated a northward flow 65% of the 
time and southward flow 35% of the time. Reise (1973), using drift-bottles over a 14 
month period in 1963 and 1964 in the southern bight, and Mcorny (1973), employing 
drogues over a two day period in August 1972, both determined surface flow within the 
bay to be generally cyclonic in nature (Figure 5 on page 10). 

Pirie and Stellar (1974) inferred surface flow within the Monterey Bay and 
vicinity through sediment distribution, obtained from visual satellite imagery and aircraft 
observations. Their flow patterns suggested large current variability but seasonal trends 
indicated northward flow along the coast from October through February’, with south¬ 
ward flow predominating from March through July. Northward flow appears to corre¬ 
late well with the classic Davidson Current regime from October through February, with 
southward flow occurring during the typical coastal upwelling period. Breaker and 
Broenkow (1989), using AVHRR imagery taken during the upwelling spring transition 
in March 1980 observed cold, upweiled water extending along the entire central 
California coast fiom a strong thermal front at Pt. Pinos. During this event, sea surface 
temperatures at Granite Canyon (located between the M onterey Peninsula and Pt. Sur, 
see Figure 1 on page 2) dropped 4°C over a six day period (Breaker and Mooers, 1986). 
Coincident with this event, current direction was observed to reverse from north to south 
in the vicinity of Pt. Sur (Wickham et al., 1987). Other satellite imagery reviewed by 
Breaker and Broenkow (1989) suggest instances of warm water intrusions entering the 
bay from the open sea. They concluded that during strong upwelling events, cold water 
was often adveacd across the entrance of the bay, then south towards Pt. Sur, forming 
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Tigure 5. Surface density analysis and drogue trajectories (21-22 Aug 1972): Note 
inferred cyclon.. circulation (Source: Moomy, 1973). 

a continuous cold frontal feature. Another common characteristic found in satellite 
images was that generally warmer waters were found inside the bay providing support 
to Bigelow and Leslie's (1930) findings that local heating plays a large role in modifying 
water masses within the bay and coastal upwclling occurs on an infrequent basis within 
the bay itself (Breaker and Broenkow. 1989). 
c. Flow at Mid-Depths 

Flow at intermediate depths (30 to 120 m) was investigated by Lammers 
(1971), using data acquired from Skogsberg, 1929-1933; Bolin ct al., 1951-1955; and the 
CalCOFI program, 1954-1967 to construct a gridded temperature field over Monterey 
Bay. Gcostrophic flow patterns wcic then implied thiough analysis of monthly-mean 
isotherm topographies (Figure 6 on page 11). '1 he study concluded that, in the mean, a 
warm core antie>clonic feature within the thcimoclinc persisted over the center of the 
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rigure 6. Monthly-mean isotherm topographies (in): Figures were obtained from 
30 ycais of data and show a highly developed thermal high situated over 
central Montcicy Bay (Source: Adapted From hammers, 1971). 


bay and tended to intensify duting the winter. This waun core feature is in contrast to 
Bigelow and Leslie's (1930) findings of slightly depressed temperatures over the canyon 
during a one week period in summer but is in agreement with hammers" inferred 
anticyclonic circulation. One possible mechanism thought to be responsible for this 
feature is momentum and heat transfer from the northward flowing CUC (Breaker and 
Broenkow, 19S9). 

d. Deep Flo Studies 

A series of deep flow studies were undertaken from 1965 to 1975 by students 
at the Naval Postgiaduatc School (NFS). These studies wcic primarily conducted near 
the canyon head at depths ranging from 91 to 485 m (T able 1 on page 12) and were 
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Table 1. SUMMARY OF NPS NEAR-BOTTOM CANYON CURRENT 


STUDIES. 


Reference 

Bottom 

Depth(s) 

Duration 

Comments 

Gatje and 
Pizinger, 1965 

130 m 

77 observa¬ 
tions on nine 
different days 

Currents tended to follow the 
canyon axis. Downcanyon flow 
was associated with incoming 
tide and upcanyon flow with 
outgoing tide. 

Dooley, 1968 

104-165m 

5-162 h 

Flows were predominantly in 
direction of canyon axis. Re¬ 
versals in flow direction showed 
a strong semidiurnal compo¬ 
nent. Sudden increases or bursts 
in current speed also were 
noted. Spectral analysis re¬ 
vealed 6 and 4 h peaks in addi¬ 
tion to semidiurnal component. 

Xjus, 1968 

150-202 m 

168 h 

Current direction was generally 
along the canyon axis. Re¬ 
versals in current direction oc¬ 
curred primarily at high and low 
tide. Spectral peaks occurred 
at 12.5, 6. and 4 h. 

Caster, 1969 

91-366 m 

3-171 h 

| 

Cross-canvon flow was fre¬ 
quently observed. Flow at an 
adjacent shelf location showed 
"eaker tidal oscillations and 
’enerallv slower speeds. 

Hollister, 1975 

485 m 

262-371 h 

Flow at 30 m above bottom was 
predominantly along canyon 
axis and apparently tidally 
driven. Strong cross-canyon 
flow was observed at 60 m 
above bottom and was poorly 
correlated with flow at 30 m. 


Source: Adapted from Breaker and Broenkow, 1989. 


primarily of short duration ranging from a day (Gatje and Pizinger, 1965) to approxi¬ 
mately two weeks (Hollister, 1975).. Most current meter depths ranged from 5 to 16 m 
above the bottom (Hollister's study used meters placed 30 and 60 m above the bottom). 
Data acquired near the canyon head revealed predominantly along-axis flow with 
semidiurnal periods. Castor (1969) and Hollister (1975), on the other hand, observed 
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frequent cross-canyon flow in the deeper, wider areas of the canyon. The spectral ana¬ 
lyses done by Dooley (1968) and Njus (1968) indicated super-tidal periods at 6 and 4 h 
in addition to the expected semidiurnal period. Gatje and Pizinger (1965) correlated 
downcanyon flow with the surface flood tide and upcanyon flow with the ebbing tide. 

Shepard (1979) conducted extensive research on flows in numerous subma¬ 
rine canyons, including the Monterey Canyon. His work in the Monterey Submarine 
Canyon was at much greater depths (up to 1445 m) than those undertaken by NPS 
(Figure 7 on page 14). Of the nine records obtained, five showed virtually no directional 
correlation with the canyon axis, unique among all of the canyons investigated. One 
suggestion for the cause of this lack of directional control was related to the meandering 
of the axis along with flows into and out of canyon tributaries. Another interesting as¬ 
pect noted by Shepard (1979) was that, unlike other California stations, the net flow was 
upcanyon in 6 out of the 9 records. Shepard stated that current direction and surface 
tides were not always coherent (Figure 8 on page 15), at times in phase with the tidal 
cycle, at other times out of phase. By observing discrepancies noted when time-velocity 
plots of different stations were overlayed (in particular stations 34 (at 357 m) and 35 (at 
384 m) and stations 58 (at 1061 m) and 59 (at 1445 m) (see Figure 7 on page 14)), then 
shifting the profiles in time to match the observed up and down-canyon flows, Shepard 
concluded that internal wave patterns were propagating upcanyon between stations at 
shallower depths but downcanyon between stations at greater depths. 
e. Internal Waves 

Almost every study done in Monterey Bay cites evidence indicative of 
internal waves. Shea and Broenkow (1982) investigated the occurrence of semidiurnal 
internal waves using hydrographic and thermistor data and resolved internal waves with 
heights of 50 to 120 m (Figure 9 on page 16). Spectral and cross-spectral analysis of 7 
days of surface tidal heights and thermistor data acquired near the canyon head sug¬ 
gested that these waves have a semidiurnal periodicity but lag the tides by about 7 h at 
the semidiurnal period. They cited the Monterey Submarine Canyon topography as the 
probable generation mechanism. As the semidiurnal long wave field approaches the 
shoaler and narrower canyon head, the period remains constant and the amplitudes in¬ 
crease as the wave energy becomes focused, as small amplitude wave theory predicts, and 
assumes the characteristics of an internal tidal bore. At high internal tides they found 
evidence that denser water from below the main thermocline within the canyon moved 
up and onto the shelves and then became pinched off as the internal tides fell 
(Figure 10 on page 17). Baroclinic tidal currents were also investigated by Broenkow 
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stations 34, 35, 58, and 59 were used to imply internal wave propagation 
(Source: Shepard, 1979). 


and McKain (1972) who found that tidal amplitudes generally decreased seaward from 
the canyon head and noted that the baroclinic and the barotropic tides were generally 
not in phase. 

2. Numerical Studies of Monterey Bay 
a. Tidal A [odds 

Schomakcr (1933) used an implicit, two-dimensional (1 km grid) program 
with realistic topography to model homogenous barotropic tidal forcing in the bay. 
Resulting height fields showed the clear progression of the tides into the bay; however, 
the inferred currents into and out of the bay associated with these heights were relatively 
weak. Consequently, other effects, such as wind forcing and ofl'shoie currents, were 
determined to dominate the flow field, rather than barotropic tidal influences. 
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Figure 8. Aloug-axis current \el«city vs. time in the Monterey Canyon: Figures 
show the incohctcnce between tidal height (smooth line) and current 
velocities (shaded). Current meter locations arc on the previous figure 
(Source: Shepard, 1979). 


b. Mean Flow Models 

Using a one-layer baiotiopic numerical cavity flow model, with boundaries 
representative of the bay, Gaieia (1971) attempted the first numerical study of the bay. 
The model was forced by momentum transfer (frictional shear) from northward flow 
adjacent (west of) to the cavity and included the effects of bottom topography 
(smoothed version of actual topography), friction, and the Coriolis force. Model out¬ 
puts were in the form of mass transport stream functions. The results of this model led 
to the conclusion that the presence of a submarine canyon can induce a closed circu¬ 
lation within the bay and that the ^-effect is negligible (Figure 11 on page 18). He also 
noted that the effects of bottom friction on the flow were negligible within the canyon 
itself but significantly influenced flow in the shallower regions. The most important re¬ 
sult obtained fiom this model was the fact that the rotational nature of flow within the 
bay was changed by suppressing or adding the bottom topography of the Monterey 











Hgure 9. Time series of temperature distributions near the canyon head: Data 
were collected 1 km (top) and 4 km (bottom) west of the canyon head 
and imply internal wave heights or 50 to 120 m. Isotherms in 'C 
(Source: Shea and Brocnkow, 1982). 


Submarine Canyon. For cavity flow withou' canyon topography, the northward mass 
tiansport stream simply entered the bay from the south and exited to the north resulting 
in cyclonically curved (low stemming from Coriolis foicc considerations. When canyon 
topography was added, momentum transfer, through shear stress from the north flow 
adjacent to the cavity, drove an anticyclonic gyre which effectively filled the Monterey 
Bay (Figure II on page 18). 

Using output from a two layer, semi-implicit, primitive equation model and 
forcing across a north-south oriented boundary, Bruner (1988). showed that the upper 
layer was insensitive to bottom topography and to lower layer boundary forcing for ve¬ 
locities less than 5 cm .s '. When lower layer flow was increased above 5 cm the 
upper layer basically mirrored the lower. The lower laver was, of course, highly sensitive 
to bottom topogiaphy. Also, when the inflow width was decreased along the northern 
boundary, and southern boundary outflow was increased, the surface flow' changed dra- 
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Figure 10. Conceptual model (top) of temperature distributions (bottom) near the 
edge of the Monterey Canyon: Top figures show the simple 2-level 
model used by Shea and Brocnkow (1982) to describe how denser water 
from within the canyon is lifted out of the canyon at high internal tide 
(a), then pinched off onto the shelves during low internal tide (b). 
Bottom figures show actual temperature cross sections along the shelf 
during high internal tides (a) and low internal tides (b) (Source: Shea 
and Brocnkow, 1982). 
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Figure 11. Garcia's one layer cavity model: Top figure represents model output 
with no canyon topography; bottom figure shows output with 
topographic efiects included. Contours are volume transport stream 
functions with velocity arrows overlaid (Source: Garcia, 1971). 

matically from cyclonic over the inner bay to totally anticyclonic flow while the lower 
layer flow was inhibited from entering the bay and turned, instead, almost immediately 
to the south of the bay where ii exited (compare Figure 12 on page 19 and Figure 13 
on page 20). Model outputs also showed that flow within the bay was highly sensitive 
to vertical shear. When no shear was present between the layers (barotropic case) the 
bottom layer had a great effect on the surface layer. However, when vertical shear was 
included (baroclinic case) the layers were essentially decoupled. Finally, Bruner noted 
that, when inflow was forced over the southern boundary simulating inflow of the 
Davidson Current, model runs depicted cyclonic circulation within the bay (Figure 14 
















Figure 12. Bruner's two-layer, baroclinic model output #1: Surface flow (top) and 
lower layer flow (bottom) with forced inflow on the northern boundary 
and outflow over the southern boundary., North is to the left. Isotachs 
arc overlaid with velocity vectors (Source: Bruner, 1988). 

on page 21). These results contradict Garcia's barotropic model outputs, which pre¬ 
dicted anticyclonic circulation within the confines of the bay from an adjacent northward 
flow. Boundary forcing in both Garcia's and Bruner's models was purely hypothetical, 
intended solely to test model responses (and the physics involved) to various flow sce¬ 
narios. 

Two final model studies, although not conducted exclusively for Monterey 
Bay, arc mentioned here because of their obvious relevance. Both models, the first 










Figure 13. Bruner's *«o layer, baroclinic model output #2: 'lhe same as in 
Figure 12 on page 19 except that inflow area was constricted and out¬ 
flow area was widened. North is to the left (Source: Bruner, 1988). 

solved analytically and the second numerically, attempt to represent the adjustment 
process that takes place when a gcostrophically balanced flow over a shelf encounters a 
submarine canyon, and arc from Klinck (1988 and 1989). 

In the first model. Klinck (1988) initially solved the shallow water equations 
by the Laplace Transform method for a two layer fluid on an /-plane, and then applied 
the boundary conditions of a flat bottom ocean intciscctcd by an infinitely long rectan¬ 
gular trench. He then applied the initial condition of an arbitrarily wide flow in the 
upper layer perpendicular to the canyon with no initial flow within the canyon itself. 









Figure 14. Bruner's two-layer, baroclinic model output #3: Surface flow (top) and 
lower layer flow (bottom) with forced inflow on the southern boundary, 
outflow over the northern boundary. North is to the left (Source: 
Bruner, 1988). 

The solution to this problem yielded standing wave motion over the canyon with waves 
radiating from cither side. In the absence of friction, these waves would continue to 
oscillate indefinitely as standing gravity modes. Results also showed that, in the steady 
state, there was no net flow across the canyon, implying that submarine canyons may 
inhibit barotropic gcostrophie flow. For shelf currents on the order of the external 
Rossby radius of deformation {R d0 ), these oscillations have periods on the order of 0.1 
to 1.0 times the local inertial ficqucncy. For "narrow" canyons (defined in this case as 












a canyon whose width is less than one half of the smaller of either the width of the cur¬ 
rent or R d0 ), Klinck showed that these oscillations would have periods on the order of 
0.2 times the inertia period. He notes that these frequency ranges are observed in most 
canyons. Breaker and Broenkow' (1989), citing this work, calculated a period of 4 h for 
Monterey Bay and noted the close agreement with Dooley (1968) and Njus' (1968) ob¬ 
served 4 h periods. They concluded that both inertial and non-linear effects associated 
with the interactions of the internal baroclinic tides and steep topography are both 
possibly at work in generating supertidal oscillations in Monterey Bay. 

Klinck (1989) conducted further experiments on the geostrophic adjustment 
process, this time utilizing a three layer numerical model (Figure 15 on page 24) where 
both the internal and external modes can act (the two top layers modeled shelf flow, the 
third layer modeled flow' within the canyon). As he states, the importance of using a 
three layer, stratified fluid is the fact that pressure gradients at the bottom of the upper 
layers will now influence the fluid within the canyon and that,".. .density redistribution 
can cancel the effect of a surface slope such that there is a steady geostrophic flow," re¬ 
sults which were not possible in his original model. In addition to the coastal current 
width, canyon width, and R the first internal Rossby radius of deformation(R dl ) now 
becomes another important controlling factor in the geostrophic adjustment process. 
Klinck initialized this model in much the same w'ay as his barotropic model, i.e. 
geostrophicallv balanced, steady, barotropic 10 cm s~ l flow over an infinitely wide flat 
bottom in the positive x direction (this would be the same as looking shoreward at a 
cross section of Monterey Bay from sea with an initial shelf current flowing equatorward 
aiong the coast) which crosses a rectangular canyon at right angles. Flow in the canyon 
was initially at rest. A 50 km wide current was used and runs were made using varying 
canyon widths including 200 and 10 km. In the 200 km wide case, the induced disturb¬ 
ance at the canyon walls are independent of each other (Figure 15 on page 24). Flow 
over the shelf decays away from the wall on the order of the current width, while flow 
within the canyon decays away from the wall on the order of R dl for the internal mode 
and on the order of the current width for the external mode. Also, the largest density 
gradients are located within the canyon which, as Klinck points out, indicates that the 
strongest currents are along the walls of the canyon and are bottom trapped. He ex¬ 
plains these results in terms of vortex stretching, stating; 
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During adjustment, the water in the canyon is pushed down the initial pressure 
gradient which causes upwelling at places where the surface elevation in the initial 
state was lowest and downwclling where it was highest. Upwelling stretches the 
vortex tubes in the canyon creating cyclonic vorticitv at the edges of the canyon. 
The mass redistribution reduces the slope of the free surface creating cyclonic 
(anticyclonic) vorticity in level 1 on the shelf side of the canyon wall in places where 
the free surface rises (falls). The isopycnal surface over the canyon near the wall 
rises sharply compressing vortex tubes creating anticyclonic vorticity right at the 
wall. 

In the narrow canyon case (10 km wide), the external mode is inhibited because the 
width is now smaller than any of the decay scales (current width, R d0 , or R dl ). Isopycnals 
within the canyon create vortex stretching (positive vorticity) resulting in a 
topographically trapped cyclone within the canyon and vortex compaction (negative 
vorticity) which results in anticyclonic rotation in the fluid above the canyon. Finally, 
Klinck noted that as the canyon width continued to be shortened (narrowed), it ceased 
to have much of a perturbing effect on shelf flow when it reached approximately one half 
of R dl . Limitations, however, should be kept in mind when applying these model results 
to Monterey Bay. First, the Monterey Bay Submarine Canyon is roughly triangular in 
shape (in plan view and in its vertical cross-section); second, the shelf region north and 
south of the bay is relatively narrow (approximately 6 km to the north and only 2 to 3 
km to the south), and; finally, the canyon is located in a bay and does not simply extend 
seaward across a coastal shelf. 

Typical values of the first internal Rossby Radius (R dl ) were calculated for 
the M ontercy Bay using the relationship: 


(*//,) 


1/2 


^cin ' j• 


where II ,, the equivalent depth, is: 


ii = *W- 

ll e 2 2 
gn n 


and X, the Brunt-Vaisala frequency, is defined as: 

v 2 _ jlIr. 

& * A A 

P 0 cz 


( 3 ) 


( 4 ) 


( 5 ) 


23 






Figure 15. Geostrophic velocity adjustment over a submarine canyon: Side view 
of Klinck's 3-levcl canyon model (top). Bottom figures show cross 
canyon (x-direction) outputs forced by a 50 km wide current over 200 
km (bottom-left) and 10 km (bottom-right) wide canyons (Source: 
Klinck, 19S9). 


















Substituting equation (4) into equation (3), and solving for the first internal mode (n= 1) 
yields: 




MI 
Kf 


( 6 ) 


Brunt-Vaisala frequencies were obtained through 2 m finite differencing techniques and 
averaged throughout the water column to obtain values for A', which ranged from ap¬ 
proximately 11 X 10~ 3 Hz over the shelf area (very stable) to 4 X 10 -3 Hz over the can¬ 
yon axis (Appendix D). Values for R d „ obtained from equation (6), ranged from 2 km 
over the shelf regions to 15 km over the axis of Monterey Bay (Appendix D). Canyon 
widths range from approximately 10 to 15 km at the mouth of the canyon, 7 to 9 km 
across the center, and 2 to 3 km near the head of the canyon. Calculated values of R 
obtained using the relationship: 




r > 
J 


( 7 ) 


ranged from 200 km over the shelf regions to 1000 km near the center of the bay (Ap¬ 
pendix D). These values are in close agreement with results obtained by Breaker and 
Broenkow (1989) through the averaging of 3 months of hydrographic data. Thus the 
Monterey Canyon could be defined as "narrow" in both the barotropic and baroclinic 
sense, although portions of the bay do approach the width scale of R dx . 

C. PURPOSE OF STUDY 

In light of recent sewage and agricultural run off problems experienced within the 
Monterey Bay, a solid grasp of the mean circulation is essential. This study attempts to 
quantify the mean flow field experienced within the bay over a four day period in May 
1988 (upwelling period) using conductivity, temperature, and depth (CTD) data, acoustic 
Doppler current profiler (ADCP) measurements, current meter data acquired by the 
Monterey Bay Aquarium, and available NOAA-11 Advanced Very High Resolution 
Radiometer (AVHRR) satellite sea surface temperature imagery. Phase relationships 
between the tides and the baroclinic flow field, along with suggested internal wave fields 
in the Monterey Submarine Canyon, are also examined. 

Finally, the circulation response within the bay to forcing by the observed mean flow 
determined from the CTD and ADCP data will be examined using an updated version 
of the baroclinic model used by Bruner (1988). Wind data, acquired both before and 
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during the period of this study, will also be used to force this model, a factor not previ¬ 
ously considered. 
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II. DATA ACQUISITION AND PROCESSING 


A. COLLECTION 

1. May 1988 NPS Student Cruise 

The vast majority of the data utilized in this study were acquired during the 
second half of an NPS student cruise conducted within the Monterey Bay during the 
period 08 to 11 May 1988 onboard the Research Vessel Point Sur, homeported at Moss 
Landing, California. Two cross-sections were run perpendicular to bottom topography, 
the first running roughly north-south across the mouth of the bay and, the second, 
roughly parallel but inshore of the first. Each line contained CTD stations which were 
sampled around the clock (Figure 16 on page 28). The seaward section intersected the 
mouth of the Monterey Submarine Canyon with station 7 situated directly over the axis 
(hereafter referred to as "Line 1"), while the landward section crossed the Monterey 
Canyon axis at station 13 and the axis of the Soquel Canyon at station 15 (hereafter 
referred to as "Line 2"). The first cross-section contained 10 CTD stations, each sampled 
a total of five times, while the shoreward Line acquired data from 7 hydrographic 
stations, each sampled four times, resulting in the overall collection of 78 CTD casts at 
17 stations throughout the period. The cruise was designed, in part, to determine the 
mean currents. Thus, to negate the effects of the tidal motions in the water column, 
each station was occupied at approximately 15.5 hour intervals, resulting in exactly a 
125% phase shift of the semidiurnal tidal cycle each time the station was sampled (Fig¬ 
ure 17 on page 29). Each station was therefore occupied on all phases of the surface 
tide. 

In addition to hydrographic data, the ADCP was operated continuously for 
each line throughout the duration of the cruise. The RV Point Sur is outfitted with an 
RD Instruments DR0150 ADCP which utilizes a four beam JANUS array operating at 
a frequency of 150 kHz. Beams are aligned fore and aft, and port and starboard, of the 
RV Point Sur. Raw data were acquired in 128 four meter bins averaged over a three 
minute ensemble interval. The ship's position was recorded from Loran-C at the end 
of each ensemble. 

2. Other Supporting Data 

Cloud free NOAA-11 AVHRR satellite imagery was obtained from 4, 8, and 11 
May 1988 and processed at the Nasal Oceanographic and Atmospheric Research Lab- 
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Figure 16. CTD stations occupied during May 1988 NPS student cruise: Stations 
1-10 (designated "Line 1") were each sampled 5 times and stations 10-17 
(designated "Line 2") were each sampled 4 times. The sampling interval 
was approximately 15.5 hours. 

oraton (NOARL) - West, located in Monterey, CA. GOES - West visual and infra-red 
imager}' were also obtained for the period of the cruise. Finally, synoptic weather charts 
were collected from the National Meteorological Center (NMC) for the months of April 
and May 1988 via facsimile broadcast. 

A continuous (5 minute averaged) record of sea surface temperature, wind speed 
and direction, tidal height, and current speed and direction was obtained from the 
Monterey Bay Aquarium (MBA) data base for the period 21 April through 20 June 1988 
(Figure 21 on page 37). The MBA weather station is situated on the roof of the 
Aquarium which is located on the shore of the southern bight at 36° 37. F N, 12F 54.0" 
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Figure 17. CTD sampling interval overlayed on observed surface tidal 
heights: Each station is represented by a separate symbol. Stations 
1*10 represent Line 1, stations 10-17, Line 2. The tide data were ob¬ 
tained from the MBA data base. 

W at an elevation of 21 m above sea level. The MBA current meter, an electromagnetic 
IntcrOcean Model S4, is situated approximately 305 ni seaward of MBA at a depth of 
6.1 in below Mean Lower Low Water (MLLW) at 36° 37.3' N, 121° 53.6' W. Tide data 
were obtained from a tide gauge located on Monterey Wharf #2 (36° 36.3' N, 121° 53.5' 
W). The gauge is operated and maintained by the National Oceanic and Atmospheric 
Administration (NOAA). Raw tidal heights were corrected to MLLW by subtracting a 
2.5 m offset. 

B. PROCESSING 

1. Neil Broun Instrument Systems (NBIS) Mark IIIB CTD Probe 

Casts were made at each station utilising a NBIS Mark IIIB CTD probe with 
manufacturer sensitivities listed in Tabic 2 on page 30. Water samples were taken at the 
bottom of each cast and sea surface temperature measurements were obtained from 
bucket measurements. Pressure offset (the pressure recorded by the CTD while sitting 








on deck) values were recorded just prior to deploying, and upon recover)’ of, the instru¬ 
ment. Finally, in order to insure accurate salinity measurements, the conductivity probe 
was rinsed with fresh water and covered after each hydrocast. 


Table 2. NBIS MARK IIIB CTD MANUFACTURER SPECIFICATIONS 


Variable 

Range 

Accuracy 

Resolution 

Pressure 

0 to 3200 dbar 

+ 3.2 dbar 

0.05 dbar 

Temperature 

-3 to +32 0 C 

+ 0.005 0 C 

.0005 0 C 

Conductivity 

1 to 65 mmho 

± 0.005 mmho 

0.001 mmho 


Calibrations were conducted prior to the cruise by checking the CTD 
conductivity, temperature, and pressure readings against standards in the laboratory. 
Differences obtained in this manner were then averaged and fit to a linear regression 
scheme in order to obtain the coefficients necessary to adjust the measurements made 
by the CTD to the reference standards. Temperature coefficients had a slope of 
0.999363 and an intercept of 0.003435 while the pressure slope was calculated to be 
1.000000. 

Post-analysis of data in the laboratory consisted of first measuring the salinity 
from the collected water samples to obtain the "true" salinity, then determining the dif¬ 
ference between salinities derived from the CTD data and those obtained in the labora¬ 
tory (Figure 18 on page 31). Points more than two standard deviations from the mean 
(outliers) were considered to be the result of random errors (human errors, ship motion 
in areas of large salinity gradients, etc.) and were removed. The mean and standard de¬ 
viation were then recalculated and outliers greater than two standard deviations from the 
mean were once more removed. A linear fit to the data obtained by this method had a 
slope of 1.036714 and an intercept of -0.347450 which resulted in a salinity offset of 
0.868 psu. Temperature and pressure values were once again compared to values ob¬ 
tained by laboratory instruments. Finally, a program which used both pre-cruise and 
post-cruise calibration coefficients, including recorded pressure offsets, was applied to 
t'nc raw' data acquired by the CTD to obtain the values used in this study. A more de¬ 
tailed description of these calibration procedures is provided by Tisch (1990). 

2. Acoustic Doppler Current Profiler (ADCP) 

After calibration of the raw data and correcting for the ship's velocity, current velocities 
to a depth of approximately 400 m were calculated. During collection, the ADCP data 
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Figure 18. Scattergram of initial salinity difference for STMAY88: Plot repres¬ 
ents difTerences between C I'D and water sample salinity at the bottom 
of each cast. Positive values indicate the instrument is reading too low. 
Mean (n) and standard deviation (sd) are indicated (Source: Tisch, 
1990). 

were averaged and stored every 3 minutes along with the ship’s position. In this way, a 
unique position is associated with each ensemble average. To recover cross sectional 
velocity profiles, data is again averaged at specified increments along the entire line. 
Inaccuracies inherent in this method include errors induced through bad navigation data 
and through the averaging techniques. For instance, if the ship was stationary at one 
of the specified increments, but underway at the next increment, the previous station 
contains more data than the latter, so averaging is somewhat biased in this respect. For 
this reason the ADCP data acquired during this study was averaged at 4 km intervals, 
which approximately corresponds to the hydrographic station spacing, to minimize both 
the errors induced through navigation, and unequal ensemble averaging. A more de¬ 
tailed description of converting raw ADCP relative velocities into meaningful absolute 
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velocities is given by Kosro (1985). Accuracies in current velocities obtained using this 
method have been cited to be on the order of 4 to 5 cm s~ l (Kosro, 1985). 

Built into the ADCP system software are constraints which determine whether 
data is to be considered "good" or not, namely: 

• The shallowest bin is no good due to side lobe interference reflected from the ship's 
hull and the ocean surface. 

• There must be at least a 60% return in signal strength. 

• Error velocity, defined as the difference in vertical velocities between orthogonal 
pairs of beams, must be within specifications in order for data within a certain bin 
to be labeled "good." 

• If the Automatic Gain Control (AGC), or reflected energy, increases greater than 
10%, the system classifies this depth as the bottom and negates data below this 
level. Typically, if no bottom is present, the AGC signal will simply decrease with 
depth (beam distance). In contrast, if a bottom is present within range of the 
ADCP beam, a jump, or increase in the AGC signal will result and data below this 
depth will be flagged as invalid. 

Data from the deepest good bin, determined by the above constraints, to the ocean 
surface is considered "good." During this cruise, ADCP data were suppressed below 
approximately 75 to 100 m in the areas near the northern half of the Monterey Canyon 
due to the AGC constraint. During replay of the acquired ADCP data over this area, 
it was noted in three out of the five passes, that a strong anomalous "spike" occurred in 
the AGC plots at the depths mentioned. As expected, the AGC plot showed gradual 
decreasing values from the surface to approximately 75 m where a sharp increase was 
noted. Below this depth, values again gradually decreased. This spike was assumed to 
be the product of biological activity of unknown origin which persisted near the bottom 
of the thermocline in these areas throughout most of the cruise. At depths below the 
AGC "spike," data exhibited low propagation loss (greater than 60% return signal 
strength) lending some credibility to the fact that data acquired below this layer may, in 
fact, have been good, although to what extent it was biased or influenced by the "scat¬ 
tering" layer is unknown. At any rate, ADCP data in these areas were re-processed by 
Mr. Paul Jessen. an NPS Oceanographer, with the AGC constraint suppressed, to re¬ 
cover the velocities in the area along the northern wall of the Monterey Canyon. 



III. METHODS AND RESULTS 


This chapter describes the results and observations obtained through the analysis 
of collected data. Methods used in obtaining stated results are provided in a separate 
subsection unless they are self-explanatory', such as in the section titled "synoptic 
meteorological situation." Evaluation of each data type is discussed separately with 
some limited comparison of results provided. An effort to tie all results together is given 
in Chapter IV. Results are presented in the following order: internal waves inferred 
through CTD data, mean geostrophic velocities, mean ADCP velocities, description of 
NOAA AVHRR sea surface temperature imagery', analysis of the synoptic 
meteorological situation both prior to and during the May 1988 student cruise, and 
finally, the methods and results of applying the mean flow field to a primitive equation 
numerical model of the Monterey Bay are presented. Some comparisons between the 
MBA data and results obtained from this study are made, when applicable. 

A. INTERNAL WAVE FIELD 

1. Methods 

To determine the extent of aliasing by internal wave processes and relate them 
to surface tidal oscillations, various depictions of the fluctuations of the IO C isotherm 
were generated. The 10'C isotherm was chosen because of its apparent proximity to the 
center of the thermocline. Data were taken from CTD measurements. As a conse¬ 
quence of the sampling interval (approximately 15 h), only waves with periods greater 
than 30 hours (frequencies smaller than the Nyquist frequency) could be resolved. The 
total structure (period, wave number, and phase speed) of the internal wave field, 
therefore, could not be determined from these data. First, a graph of the 10‘C isotherm 
pressure across the mouth of the bay (Line 1), along with the mean isotherm pressure 
was generated (Figure 19 on page 34). To obtain a relationship between surface tidal 
height and this internal wave field, plots of the surface tidal height versus isotherm 
pressure were then produced for each station (Figure 20 on page 35). 

2. Results 

The plots of the 10'C isotherm pressure changes during each pass of Line 1 
(Figure 19 on page 34) show the active internal wave structure encountered during data 
collection. The largest amplitudes (up to 30 m) were observed to be in the vicinity of the 
shelf breaks (stations 3 and 9 at 9 and 30 km on flic graph), over the axis of the 
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Figure 19. Isotherm variance along Line 1: Top figure shows the pressure of the 
10’C isotherm across the mouth of the Monterey Bay during each pass 
(see legend). North is to the left, south is to the right. Bottom figure 
shows the mean 10”C isotherm depth and the calculated isotherm vari¬ 
ance at each station along the line. Note the large amplitude fluctu¬ 
ations near the shelf break (variances of 94.36 m 2 and 69.43 m 2 ) and 
over the center of the Monterey Submarine Canyon (variance of 78.14 


m 2 ). 
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change in pressure of the IO'C isotherm as a function of the observed 
surface tides, ’fide data were obtained from MBA data base. 


Monterey Submatine Canyon (near the 24 km mark on the graph), and near shore over 
the shelf at station 1. What appear to be nodal points can be seen near stations 2 and 
6. The plotted mean 10‘C isotherm pressuic (Figure 19 on page 34) shows gradual 
deepening of the isotherm from the northern part of the bay (station 1) to the center of 
the canyon (station 7) then steep shoaling from station 7 to the southern shelf at station 
10. Calculated \ariance shows high values in the aieas of large amplitude fluctuations 
picviously mentioned. In general, the observed internal wave amplitudes were smaller 
in Line 2 (stations 11-17) with most amplitudes falling within the 10 to 15 m range. This 
was somewhat surpiising considering past studies (Shea and Broenkow, 19S2, and 
Biocnhow and McKain, 1 ( )72) which showed highet amplitude internal waves nearer the 
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head of the canyon than towards the mouth. This again may simply be the result of the 
relatively large station sampling interval, i.e. the internal waves in Line 1 may have had 
a different frequency or phase relation with the surface tidal heights than those en¬ 
countered in Line 1, resulting in poorer resolution of these waves in Line 2 data. 

The 10’C isotherm pressure versus tidal height at each station of Line 1 
(Figure 20 on page 35) show no direct surface tidal correlation. Large differences in 
pressure can be seen even though the surface tidal height is at nearly the same level. 
For instance, at station 8, when the tide is approximately 0.5 m, the pressure of the 
10'C isotherm ranges from 20 to 45 dbar (note 1 dbars:l m in these ranges). Similar 
cases of varying isotherm pressure at constant tidal heights can be seen at other stations 
(station 1 had isotherm pressures which varied by as much as 25 dbar during two sepa¬ 
rate periods when the tide was at the same height). If internal tides were in phase with 
the external tides, these graphs should basically show a positive sloping line from the 
bottom left to the upper right (i.e. the isotherm would rise and fall with the external tidal 
height). Station 9, near the southern shelf break, is the only station which comes close 
to this description, where all points nearly lie on a straight line. A cross-spectral analysis 
between one week of continuous thermistor data near the canyon head and surface tidal 
heights, conducted by Shea and Broenkow (1982), revealed a 7 h phase lag. A similar 
analysis of approximately 3 months of thermistor data, from mid April to mid June 1988, 
collected from the MBA sensor at a depth of 6 m in the southern bight (Figure 21 on 
page 37), revealed an 8 h phase lag at the predominant semi-diurnal 12.8 h period. 
However, phase relationships between surface tidal heights and internal waves near the 
canyon head and over the shelf regions more than likely differ from those in the deeper 
regions of the bay, especially over the canyon. 

In summary, the internal wave field encountered during the 08 to 11 May 1988 
data collection period was highly energetic, especially along Line 1. Amplitudes varied 
across the bay, but appeared to reach a maximum over the shelf breaks which lends 
support to Shea and Brocnkow's (1982) proposed model of processes occurring near the 
edge of the Monterey Canyon (see Figure 10 on page 17). Although the direction of 
propagation of these internal waves could not be determined from the available mass 
field data, the high variability across the mouth of the bay would seem to suggest at least 
some cross canyon component exists. Possibly the most useful information gained from 
the above analysis is the fact that the isotherm slopes created by the internal waves ex¬ 
ceed the mean slopes due to geostrophic adjustment by two orders of magnitude. This 
implies that geostrophic velocities derived from one pass alone is virtually useless be- 
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Figure 21. Acquired MBA ocean and atmospheric data 21 April - 20 June 
1988: MBA data is averaged and stored ever}’ 5 minutes. Above 
curves were generated using a 3 hour running mean filter. Dashed lines 
represent total mean. 
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cause the internal wave "noise" will quite effectively mask the mean signal. Therefore, 
some sort of averaging scheme is a necessity. The next sections describe the geostrophic 
flow obtained from the mean mass field and the mean ADCP field. 

B. MEAN TEMPERATURE, SALINITY, DENSITY, AND GEOSTROPHIC FLOW 
1. Method 

Analysis of the hydrographic data acquired during the May 1988 NPS student 
cruise was conducted in two steps. First, data from each pass were run through a com¬ 
puter program which utilizes horizontal density gradients between stations via the ther¬ 
mal wind equation: 


8v 8 cp 

dz JPq cx 

to calculate geostrophic velocities throughout the water column. The second part con¬ 
sisted of averaging the mass field. This was accomplished by averaging the temperatures 
and salinities at each station (from 5 samples in Line 1 and 4 samples in Line 2). From 
this mean field, geostrophic velocities and cross-sections were then produced (using 
equation (8)) for both Lines (Figure 24 on page 42, Figure 26 on page 47, and Appen¬ 
dix A). Because of the sharp topography, the assumed level of no motion used in all 
geostrophic velocity calculations was the deepest common depth between stations. 

As previously stated, the sampling interval was designed so that tidal oscillations 
could be averaged out of the water column. To measure the success of this particular 
aspect of the cruise, the mean tidal height for the period was calculated using MBA data, 
and the variance of tidal heights (from observed heights at each station during sequential 
passes) about this mean were then derived (Figure 22 on page 39). The calculated var¬ 
iance was relatively small at 3.45 x KHm 2 for Line 1, and 1.12 x KHm 2 for Line 2 
(Figure 22 on page 39). Initially, averaging schemes using various pass combinations 
were tested in an attempt to minimize the variance of station tidal height about the 
MBA tidal mean. Although averaging all line passes did not produce the minimum 
variance about the MBA mean, it was decided to use all five because it was apparent 
that the internal tides were not in phase with the surface tidal heights (discussed in sec¬ 
tion on internal waves). 

To further measure the "success" of this averaging process, the net mean mass 
transport across each line was calculated. In the mean, the net geostrophic mass trans¬ 
port into and out of the bay should be zero, assuming that the barotropic transport was 
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Figure 22. Effective tides: 'I lie plot shows surface tidal fluctuations from 5 min¬ 
ute averages of observed tidal height obtained from the MBA data 
base. Solid line shews mean MBA tidal height, dotted line shows cf- 
fccthc surface tide derived thiough averaging of tidal heights observed 
at each hydrographic station (sec also Figure 17 on page 29). Variance 
wui. respect to the MBA mean shown at the bottom right. 


minimal, Calculated net volume tiansport was .152 Sv for Line 1 and .012 Sv for Line 
2, both considered balanced within the resolution of the methods used within this study 
for reasons discussed shortly. 

Finally, the standard error of the mean gcostrophic velocity (aj) was calculated 
(Figure 24 on page 42 and Figure 26 on page 47) using the relationship: 


°x = 



( 9 ) 


where the standard deviation (a) is: 


a — 


Z(X-Xf 

n 



( 10 ) 
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The derived geostrophic velocities (X) between each station on each individual pass was 
used as the population (n). 

In addition to geostrophic cross-sections, maps of the temperature fields at 0, 
50, and 200 dbar, dynamic heights at 0 dbar relative to 50 dbar and 100 dbar and 200 
dbar relative to 500 dbar, as well as cross sections of temperature, salinity, and density, 
were generated (Figure 27 on page 48, Figure 28 on page 49, Figure 23 on page 41, and 
Figure 25 on page 46), The dynamic height topographies were then used to infer the 
geostrophic flow patterns. 

In summary, three parameters were calculated for each data set to determine the 
effectiveness of the averaging process in extracting the mean flow:, station tidal variance 
about the MBA mean tidal height, section mean volume transport, and the standard 
error. Experimentation with the averaging of various pass combinations, to minimize 
the tidal variance of the data set about the MBA mean tidal height, did not reduce the 
net volume transport (i.e. increase the implied accuracy of the mean) through the sec¬ 
tion, This further illustrates the fact that the internal wave field was not in phase with 
the surface tidal heights. If they were in phase, the most accurate mean flow (smallest 
net volume transport) would naturally correlate with the averaging scheme which yielded 
the smallest variance about the mean MBA tidal height during the period. Since this 
was not the case, the decision to average all 5 passes in Line 1 and all 4 passes in Line 
2 was arrived at, from a purely statistical viewpoint, to maximize the number of degrees 
of freedom. Statistically, this method worked well since calculated standard errors were 
smaller than the derived mean. With this in mind, the overall data acquisition method 
was successful at minimizing the variance of the tidal heights and creating a nearly bal¬ 
anced mean volume transport by virtue of pass averaging alone. The method was not 
entirely successful however, at systematically (deterministically) reducing the high fre¬ 
quency internal motions. 

2. Results 

a. Mean Property Cross Sections, Geostrophic Flow, and Mass Transport 

(1) Line 1, The Line 1 mean temperature and salinity cross sections 
(Figure 23 on page 41) indicate an area of warm, fresh water at the surface (greater than 
12'C) extending across the mouth of the bay from station 5 to station 9. This warm, 
fresh "pool" of water appears well imbeded in the vertical structure and extends down 
to approximately 100 dbar, as the 8.5 and 9 'C isotherms indicate. At deeper depths, over 
the Monterey Submarine Canyon, upwelling of cold water, most notably along the 
northern shelf break, appears to dominate the temperature field as shown by the 8 and 
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Figure 23. Line 1 mean temperature, salinity, and density cross sections: Figures 
depict mean temperature (°C, top), salinity (psu, middle), and density ( 
kg nr\ bottom) sections for Line 1. Caution should be noted in inter¬ 
preting deep isolines adjacent to the canyon walls. These steep slopes 
more than likely do not exist and are merely the product of the con¬ 
touring program used to develop these graphs. Station positions are 
indicated across the top of the salinity section (sec also Figure 16 on 
page 28 for locations). 
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Figure 24. Line 1 mean grostrophic flow field and standard error: Top figure is the 
mean geostrophic flow, bottom is the standard error.. Hydrographic 
stations are notated at the top of the section and are situated across the 
mouth of Monterey Bay (Line 1, see also Figure 16 on page 28 for 
station locations). Solid lines depict flow into the bay, dashed lines 
show flow out of the bay. Geostrophic velocities were calculated using 
the deepest common depth between stations as the level of no motion. 
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8.5*C isotherms. The salinity cross section shows the fresher water (33.65 psu) located 
between stations 6 and 9. Again, some indication of upwelling of denser, saltier water 
from the canyon onto the shelf region is indicated (again most notably over the northern 
shelf) by the 34.05 and 34.00 isohalines. Density anomaly cross sections (Figure 23) 
reflect the temperature and salinity fields by indicating relatively light water located in 
the top 100 dbar of the water column and centered nearly over the canyon axis with 
denser water extending up onto the northern shelf along the canyon wall from a depth 
of 200 m. 

The calculated mean geostrophic flow (Figure 24 on page 42, and 
Appendix A) reveals two apparently separate flow patterns. At the surface (0 to ap¬ 
proximately 75 dbar) flow is split into two regimes: broad inflow across the northern 
half of the Monterey Bay from stations 1 to 7; and, outflow across the southern half of 
the bay from stations 7 to 10. Inflow extends to a somewhat greater pressure (125 dbar) 
than the outflow region (50 dbar). Maximum surface inflow was calculated at 15.6 cm 
s~' between stations 3 and 4 and a maximum outflow of -18.6 cm s-’, situated between 
stations 9 and 10, was indicated. Corresponding standard errors of the mean at these 
locations were 9.9 and 4.8 cm .s -1 respectively (Figure 24 on page 42). Analysis of the 
current speed and direction from the MBA current meter in the southern bight area 
(closest to station 10) over the period showed northwest/southeast flow reversals with 
the seaward flow predominating (Figure 21 on page 37). Also evident was the corre¬ 
lation between the strongest flow (highest velocities) and the seaward direction. Mean 
current speed and direction, derived from these data during the period of the May cruise, 
was 4.5 cm s~‘ towards 245T, or seaward. This outflow is smaller than the noted 
geostrophic flow between stations 9 and 10, not unusual considering the current meter's 
close proximity to shore. The second flow pattern is located over the Monterey Sub¬ 
marine Canyon between 100 and 600 dbar and consists of outflow along the northern 
wall of the canyon, with inflow indicated over the southern half of the canyon. The 
cores are found at 250 dbar (downcanyon flow of 10.6 cm .r 1 ) between stations 5 and 
6, and at 150 dbar (upcanyon flow of 10.1 cm s* 1 ) under station 8. Standard errors cal¬ 
culated at these positions and depths yielded values of 8.3 and 9.1 cm 5' 1 respectively, 
showing that the deep flows were barely resolved with this averaging technique. The 
mean flow pattern shows strong vertical shears near stations 6 and 8, along with weaker 
horizontal shear zones over the center of the bay. Standard errors at each level 
(Figure 24 on page 42) indicate that the derived surface mean flow is much more re- 





solved (especially the outflow in the southern half of the bay, past the Monterey Penin¬ 
sula) than the mean flow at depth, which barely surpasses the standard error values. 

The net mass transport calculated for the mean flow was .152 Sv. 
This result was not surprising and appears to be within the resolution of the geostrophic 
methods exploited for a number of reasons, namely: 

• The geostrophic method does not take into consideration any non-linear effects, 
which can produce errors oip to 10% in some areas of the ocean such as near 
eddies, high shear zones, etc. Because of the Monterey Bay's geometry, and ap¬ 
parently active internal wave system, non-linear processes more than likely are at 
work and are neglected in geostrophv. Rough estimates of the extent of non¬ 
linearity were obtained through calculations of the Rossby number (R„): 

R 0 = j[> 00 

where f the coriolis parameter is 0(10^*); velocities (V) are 0(.l m 5"‘); and typical 
length scales (L), estimated by the computed first internal Rossby radius of defor¬ 
mation (Appendix D), of 0(5 to 15 km). Values range from approximately .1 to 
.2, implying that up to 20% of flow may be affected by non-linear processes in 
Monterey Bay. 

• Recent studies (Tisch, 1990) indicate that geostrophic velocities calculated between 
stations with close spacing typically have higher errors than those calculated at 
larger grid spacing. When station spacing is less than approximately 4 to 5 km, 
which was the case in this study for nearly all the stations, small errors in posi¬ 
tioning and measurement, along with the ship drift, tend to propagate larger errors 
in geostrophic velocity calculations than values derived from stations situated far¬ 
ther apart. 

• Geostrophic velocities were calculated using the common deepest depth between 
stations. Thus, no deep flow (greater than 1000 dbar) is resolved in this study. 

• Approximately 4 km of the bay, on b ‘h ends of Line 1, were not covered. This 
was particularly important in the southern part of the bay where strong outflow 
was indicated right up to station 10. If strong outflow continued right up to the 
Monterey Peninsula, a mass transport rate of roughly -.036 Sv ( 
4000 m x 50 m x 0.18 m.r 1 ) was missed in Line 1. 

• The barotropic flow was not resolved through geostrophv. Rough estimates of this 
component were obtained by first calculating the surface area within the bay and 
bounded by Line 1 (882 km J ), then computing the flow that would be required to 
raise the water level 1.5 m (approximate tidal range of the bay) and the ensuing 
transport rate. From these calculations, barotropic flow was determined to be on 
the order of .1 to .3 cm y' winch would account for a mass transport of ± ,06 Sv. 

Therefore, errors on the order of .1 Sv can be accounted for by reasoning of missed 
coverage across the bay, deep flow, and barotropic transport. Errors of O(.l) Sv are also 
probably a realistic estimate of °rrors propagated through geostrophv (i.c. non-linearity, 
station -pacing, assumed level of no motion, etc.;. In view of this, and calculations of 
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the standard error, the derived mean flow field appears balanced within the limitations 
of available data and methods. Within the confines of this study, the term "mean flow" 
refers to the flow during the period from 08-11 May 1988 and may not represent a 
"typical" flow pattern in the Monterey Bay (if such a pattern even exists). An attempt 
will be made later to relate this "mean flow" to processes believed to have been occurring 
at the time these particular data were collected. 

(2) Line 2. The mean Line 2 temperature and salinity cross sections 
(Figure 25 on page 46) show smaller perturbations in the thermocline than the Line 1 
data. Somewhat warmer water can be seen at the surface over the northern and south¬ 
ern shelf regions but overall isotherms are relatively level. This is in contrast to the Line 
1 isotherms which show a temperature maximum over the center of the canyon. This 
may be a consequence of increased heating of the shallower shelf waters which rapidly 
dominate the bathymetry as the canyon head is approached. This warmer water situated 
over the shelf regions of the inner portions of the bay can also be seen on the OS May 
1988 satellite imagery (Figure 32 on page 58), discussed later. The salinity field also 
displays relatively level isohaline contours within the thermocline except for somewhat 
fresher water indicated over the southern shelf region (Figure 25 on page 46). The cross 
section of mean density shows somewhat lighter (less dense) water situated over the 
northern and southern shelf regions as indicated by the 25.6 isopycnal. 

The mean currents derived for Line 2 shows broad inflow over the 
center of the section between stations 15 and 12 with a maximum at the surface between 
stations 12 and 13 of 7.4 cm s~' and a subsurface maximum at 50 dbar of 6.2 cm s~' 
(Figure 26 on page 47). Standard errors were calculated at 2.7 and 1.4 cm s -1 respec¬ 
tively.. The calculated mean outflow was again narrower, but stronger, than the inflow, 
and confined mainly to the southern portion of the bay, over the shelf region, between 
stations 11 and 10. Maximum outflow was calculated at 15.9 cm s' 1 in this area with a 
standard error of 6.2 cm .r 1 As noted earlier, both mean velocities (hence, mass trans¬ 
port) and velocities calculated during each pass, were smaller than those calculated for 
Line 1, implying that much of the circulation derived in Line 1 is closed, and does not 
extend significantly into the heart of the bay itself. The net mass transport, derived from 
mean flow values, was .012 Sv (an order of magnitude smaller than that of Line 1). 
Again, similar arguments for not attaining mass balance closure in Line 1 can also be 
applied to Line 2. 

(1/ Mean Temperature and Dynamic Height Maps. The mean surface 
temperature map (Figure 27 on page 48), interpreted by hand, depicts two pools of 
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Figure 25. Line 2 mean temperature, salinity, and density cross sections: Figures 
depict mean temperature (°C, top), salinity (psu, middle), and density 
(kg nr 3 , bottom) sections extending roughly north-south across the 
center of the Monterey Bay (Line 2) and shows both the Soquel (left) 
and Monterey (right) Submarine Canyons. Caution should be noted 
in interpreting aecp isolines adjacent to the canyon walls. These steep 
slopes more than likely do not exist and are merely the product of the 
contouring program used to develop these graphs. Station positions 
arc indicated acioss the top of the temperature section (see also 
Figure 16 on page 28 for locations). 
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Figure 26. Line 2 mean geostrophic flow field and standard error: Top figure is 
mean geostrophic flow, bottom is standard error. Hydrographic 
stations are notated at the top of the section and run roughly north- 
south across the center of Monterey Bay (Line 2, see Figure 16 on page 
28 for station location). Solid lines depict flow into the bay, dashed 
lines show flow out of the bay. Geostrophic velocities were calculated 
using the common deepest depth between stations as the level of no 
motion. 
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Figure 27. Mean temperature fields: Mean temperature fields ('Q on the 0 dbar 
(top-left), 50 dbar (top-right), and 200 dbar (bottom) pressure surfaces. 












Figure 28. Mean dynamic height fields: Mean dynamic heights at 0 dbar relative 
to 50 dbar (top-left) and to 100 dbar (top-right) and at 200 dbar relative 
to 500 dbar ('bottom'). 
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warm water at the surface. The first is centered over the mouth of the canyon with a 
12.45'C maximum at station 7, and the second appears to be contained in the northern 
bight area with a 12.43'C maximum at station 17. The pool of warm water centered at 
station 7 maintains its identity down to 50 dbar (the approximate base of the mixed 
layer), with station 7 retaining the maximum 50 dbar temperature at 9.61’C. Evidence 
of the warmer water present at the surface in the vicinity of the northern bight disap¬ 
pears at depth. This may be a consequence of the previously mentioned surface heating 
of shelf waters in the northern bight, which typically does not extend down through the 
thermocline. Other satellite SST data (Tracy, 1990) also show enhanced surface heating 
in the northern bight, possibly since the area is sheltered from the northwesterly winds 
by the Santa Cruz mountains. Temperatures at 200 dbar are slightly colder over the 
canyon, especially towards the northern wall, as evidenced by the minimum temperature 
(7.65’C) at station 6. 

The mean dynamic topographies (Figure 28 on page 49) at the sur¬ 
face relative to 50 and 100 dbar, show higher relative heights near the mouth of the bay, 
roughly centered over the canyon axis with maximum values of .192 and .194 dy m at 
stations 6 and 7 respectfully and implies anticvclonic geostrophic circulation within the 
mixed layer over the Monterey Bay. Deeper circulation is cyclonic with low dynamic 
heights over the canyon axis at stations 6 and 7 of .388 and .386 dy m respectively, 

C. MEAN FLOW FROM ADCP DATA 
1. Method 

As an alternative method in determining a mean flow (and to compare with 
gcostrophically derived mean velocities), ADCP data were analyzed for Line 1. This was 
accomplished by first breaking the ADCP data up into two portions, the first with the 
ship on a heading of 180’ T (i.e. between stations 1 and 7), and the second part, when 
the ship was on a heading of roughly 145T (between stations 7 and 10). This was nec¬ 
essary so the applied program could compute velocities that would correspond to the 
directions obtained from geostrophic calculations. Data were averaged every 4 km along 
each portion of Line 1 (stations 1-7 and 7-10). After the velocities for each pass were 
computed, they were again averaged, over all 5 passes. Horizontal distance averaging 
was done since the 3 minute profiles were not co-!ocatcd from pass to pass (see 
Figure 29 on page 51). Finally, mean files from both portions (stations 1-7 and 7-10) 
were combined to produce one plot of ADCP mean velocities that would extend from 
station 1 to 10 so that direct compaiisons could be made with the geostrophic velocities. 
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Figure 29. Line 1 ADCP averaging scheme: Top figure shows averaging interval 
for stations 1-7 and the bottom figure for stations 7-10. Horizontal 
averaging as well as vertical averaging of each pass was done to mini¬ 
mize errors induced through slight offsets of data profiles along each 
pass. 
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(see Figure 30 on page 53, and Appendix B for numerical outputs). Standard errors 
were calculated using the mean ADCP velocities and the velocities obtained during each 
separate pass in the same manner as was done with the geostrophic data. 

2. Results 

a. ADCP Velocities Obtained During Successive Passes of Line l 

Contrary to velocities obtained through geostrophy from individual passes, 
ADCP velocities will not be influenced by mass fluctuations and detect the total veloci¬ 
ties (both the geostrophic and ageostrophic components) present at the time. Although 
ADCP velocities collected during each pass undoubtably are also influenced by internal 
wave motions, it is felt that results obtained during each pass warrant a brief description. 

ADCP velocities obtained during successive passes of Line 1 ranged from 
light flow up to 20 cm s _1 in pass 3. Numerous current reversals between successive 
passes were noted and effectively demonstrate the active dynamical processes at work 
within the bay. Also, large shear zones, both in the vertical and horizontal, were in ev¬ 
idence in virtually all passes. In general, ADCP velocities were less than the geostrophic 
velocities calculated during each pass, especially at depth. This is not surprising since 
the calculated geostrophic response to the large swings in isopycnal slopes, created by 
the large amplitude internal waves, would tend to produce unrealistically large velocities, 
especially in the thermocline. 

No definite relationship between the surface tidal heights and flow within 
the bay could be determined from the ADCP passes. During pass 3 of Line 1, the sur¬ 
face tide was in the ebbing phase between stations 4 and 9 (see Figure 17 on page 29), 
going from Higher High Water (HHW) at station 4 to Higher Low Water (HLW) at 
station 9. During this period the ADCP data depicted virtually all flow going upcanyon 
(into the bay) with a 20 cm maximum centered over the canyon axis at a pressure of 
approximately 100 dbar. On the next pass (pass 4), the tide was in a flooding stage be¬ 
tween station 5 (at Lower Low Water (LLW)) and station 10 (at Higher High Water 
(HHW)). In this case a total reversal in flow was noted with no evidence of the 20 cm 
y‘ maximum previously noted, Virtually all flow was directed out of the bay with a -18 
cm s'' maximum at the surface near stations 8 and 9. These two passes tend to support 
Gatje and Pizingcr's (1965) current meter observations which correlated downcanyon 
flow with surface flood tides, and upcanyon flow with surface tidal heights in the ebbing 
stage. The apparent correlation between surface tidal heights and bay flooding flushing 
periods implied by these two passes does not hold in the other passes however. For in¬ 
stance, pass 1, again in a flooding stage from LLW at station 1 to Lower High Water 
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Figure 30. Line 1 mean ADCF velocities and standard error: Top figure shows 
mean ADCP velocities along Line 1. Negative values indicate flow out 
of the bay, positive values indicate flow into the bay. Bottom figure 
represents standard error values calculated for the ADCP velocities. 
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(LHW) at station 10, depicted outflow over the northern shelf and along the north can¬ 
yon wall, but inflow over the :enter of the canyon and over the southern shelf. Also, 
during pass 2, which went from HHW at station 4 to nearly LLW at station 10 (ebbing 
tides), the ADCP velocity section depicted outflow at the surface but inflow at 100 dbar. 
These passes tend to support observations made by Shepard, et al. (1979), and by 
Broenkow and Mckain (1972), that the baroclinic and barotropic tides are generally not 
in phase. 

b. Mean ADCP Velocities and Standard Errors 

The mean ADCP velocity structure shows broad surface inflow over the 
northern half of the bay from station 1 to station 7 with the core situated approximately 
17 km south of station 1 having a maximum velocity of 15 cm s~'. A narrower region 
of outflow is revealed along the southern half of the bay, south of station 7, with a 
maximum velocity of 10.5 cm s~' centered approximately 34.5 km south of station 1 
along the Line 1 track. Weak outflow (-2 cm s' 1 ) is indicated near the northern canyon 
wall at 182 dbar with stronger inflow (13.7 cm s~' at 150 dbar) indicated over the 
southern half of the canyon axis between stations 7 and 8. Standard errors ranged from 
1.7 cm s _l for the surface inflow center to 4.1 cm s-‘ for the core of the deeper inflow 
between stations 7 and 8 (see Figure 30 on page 53). This represents a better signal to 
noise ratio for the ADCP data than for the velocities derived through geostrophy and 
implies that ADCP techniques may have an advantage over geostrophy when attempting 
to determine mean flow patterns within the Monterey Bay. 

c. Comparison of Mean ADCP and Geostrophic Velocities 

Table 3 on page 55 summarizes the maximum velocity values and their lo¬ 
cations, obtained by using both mass averaged geostrophy and by averaging the ADCP 
velocities. Considering the major differences in the two techniques, the results of both 
methods compare favorably in the .. ation and in the magnitude of the flow centers. 
Horizontal differences in the location of maximum velocities were at most 4 km's apart. 
This is expected since the ADCP continuously collects and averages velocity data, while 
geostrophic velocities can only be calculated for the water column centered between 
CTD stations. Maximum velocities at depths were also fairly close, differing by an av¬ 
erage of 50 m. Velocity magnitudes were comparable for the surface flows, differing by 
an average of 4 cm s '. At depth, velocities at the inflow region agreed well on the 
southern side of the canyon, with differences of only 3.6 cm s~' between the two meth¬ 
ods. The strong geostrophic outflow region at depth along the northern slope of the 
canyon did not fare as well. The sense of flow was the same but the velocities differed 
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by almost 8 cm srK A few possible explanations for this discrepancy in flow at this lo¬ 
cation follow: 

• Anomalously high AGC values were observed during the processing of the ADCP 
data (see Chapter II) over this region. The "scattering layer", cited as the possible 
culprit for the unusual signal, could have adversely effected ADCP measurements 
below, rendering data collected in this area "suspect." 

• Standard error calculations are smaller than mean values of both geostrophic and 
ADCP flow, with the exception of the ADCP mean obtained for this deep outflow 
region where the mean is -2.10 cm s' 1 and the standard error is 3.18 cm s _ >. This 
is another reason to label the ADCP data obtained in this area as "suspect". 

• CTD station spacing was very tight (approximately 2 km vice 4 for the other 
stations) over the northern canyon shelf whereas ADCP velocities were averaged 
every 7 4 km over the whole of Line 1. This larger averaging interval could have re¬ 
sulted in smoothing the center of the outflow core and resulted in a much smaller 
velocity. 


Table 3. SUMMARY OF LINE I MEAN VELOCITIES: Negative velocities in¬ 


dicate flow out of the bay, positive values denote flow into the bay. 


Type 

Location (km 
from sta 1) 

Depth(dbar) 


Std Err 
(cm .r 1 ) 

GHOSTROPIIIC 

13.16 

0 

15.62 

9.95 

ADCP 

17.32 

0 

14.97 

1.72 

GEOSTROPHIC 

33.17 

0 

i 

GO 

4.81 

ADCP 

34.48 

0 

-10.48 

3.49 

GEOSTROPHIC 

16.78 

250 

-10.65 

8.3 

ADCP 

17.32 

182 

-2.10 

3.18 

GEOSTROPHIC 

25.89 

150 

10.15 

9.11 

ADCP 

29.99 

118 

13.71 

4.10 


D. DESCRIPTION OF NOAA AVHRR SATELLITE IMAGERY 

Cloud-free Satellite imagery' was available for 04, 08, and 11 May 1988 (Figure 31 
on page 57 and Figure 32 on page 58). The 04 May image reflects strong coastal up- 
welling near Pt. Ano Nuevo, 35 km to the north of the Monterey Bay. This cold, up- 
welled water (7 to 8’C) appears to extend into Monterey Bay from the north. Warm 
water from inside the bay appears to be advecting southward past Pt. Sur. The implied 
circulation is inflow at the north of the bay and outflow at the south. 
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The images taken on 08 and 11 May are nearly identical and only the 08 May image 
is shown. Skin temperatures over the bay showed a sharp rise in temperature of ap¬ 
proximately 5'C over the four day period from 04 May. Whether this is a result of warm 
water advected into the bay, or simply a product of solar heating which ensued upon the 
cessation of winds and upwelling, cannot be determined. Coastal upwelling is still in 
evidence to the north and south of the bay on the 08 May image. Colder waters, either 
being upwelled along the coast or advected into the northern part of the bay, again ap¬ 
pear on the 11 May imagery. 

E. METEOROLOGICAL SYNOPSIS 

1. The Period Prior to Data Acquisition (25 April - 08 May 1988) 

On 25 April, the central California coast was dominated by a 1026 mb high 
pressure center situated over the center of the state. This pattern resulted in a relatively 
light pressure gradient overall, allowing the mesoscale land-sea breeze processes to 
dominate the flow pattern. A low pressure system was located approximately 900 n mi 
west of Washington state. By the 26 rt the high began to break down under the influence 
of this low pressure system and by the 27 rt , this low, with its associated cold front, passed 
over the Monterey Bay. While no significant weather or winds were produced by this 
system, it did clear the path for a second cold frontal system which was approximately 
600 n mi west of Monterey Bay at this time. By 28 April, the pressure field in the vicinity 
of the Monterey Bay began feeling the influence of this system generating a light (3 to 
5 m s~‘) southwesterly pre-frontal wind region regime. This strong frontal system passed 
over the Monterey Bay early in the morning of 29 April causing the winds to rapidly veer 
to the northwest and increase in magnitude to 18 m .rC As this front continued east¬ 
ward across the state, strong high pressure ridging behind this front moved over the 
Monterey Bay and maintained a strong pressure which continued generating 13 to 18 
m s~' winds over the area through the morning of 01 May. From the afternoon of 01 
May through 07 May, the winds remained out of the northwest at moderate velocities 
(5-8 m s’ 1 ) increasing slightly the morning of the 4* when a third front, weaker than the 
previous one, passed over the Monterey Bay. In summary, this period was very "unset¬ 
tled" for May, with three separate frontal systems passing over the bay. Winds were 
predominantly out of the northwest to north, at times very strong, for a period of ap¬ 
proximately nine days, from 29 April to 07 May. 
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Figure 31. 


NOAA AVHRR satellite imagery from 2336 Z 04 May 1988: Darker 
shades along the coast to the north and south of the Monterey Bay 
depict colder water. Lighter shades over the extreme eastern portion 
of the bay and within the northern and southern bights of the bay de¬ 
pict wanner water. The light streak extending southwest from the 
Monterey Peninsula represents warmer waters. The large areas of black 
over the ocean are cloud cover. Latitude and longitude lines are in 
one-half degree increments. Rectangle at left is SST scale meant for 
use with color copy and is not usable with this black and white re¬ 
production. 


57 






Figure 32. NOAA AVHRR satellite imagery from 2247 Z 08 May 1988: Darker 
shades within the Monterey Bay and to the west of the bay depict 
warmer water. Darker areas along the coast, south of the Monterey 
Peninsula, depict cooler upwelled water. Black areas over the ocean are 
clouds. Latitude and longitude lines are the same as Figure 31. Rec¬ 
tangle at bottom is SST scale meant for use with color copy and is not 
usable with this biack and white reproduction. 
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2. Synoptic Situation During the Period 08 - 11 May 1988 

The synoptic meteorological situation during the period of data acquisition was 
high pressure, centered to the southwest of the bay, which was strong enough during this 
period to block the advance of any atmospheric low pressure systems or cold fronts. 
The pressure gradient was not strong enough to influence the mesoscale land-sea breezes 
which dominated this entire period (see Figure 21 on page 37). Calm conditions pre¬ 
vailed. with very low (less than 1 m) sea states during the entire survey cruise. 

3. The Monterey Bay Aquarium Wind Field 

Review of the MBA wind data taken over a 3 month period, from the last week 
of April to the last week in June 1988, shows the predominant diurnal land-sea breeze 
fluctuations evident during periods of light winds. The wind events previously men¬ 
tioned during the last part of April and continuing to approximately 07 May is evident 
in this time series (Figure 21 on page 37). During this period winds were predominantly 
out of the northwest and attained velocities of up to 15 m 5~‘. The overall mean wind 
speed for the three month period was approximately 3.5 m ,r> from 2607'. Calculated 
mean wind speed and direction for the period of the May 19SS cruise was 2.3 m from 
260T . Also of note is the predominantly stronger winds during afternoon seabreeze 
events, with winds dying off rapidly during the evening hours. 

F. SUMMARY OF MEAN CURRENT FLOW 

Surface How implied through satellite imagery, geostrophic flow, ADCP flow, and 
dynamic height analyses were consistent, They show mean surface inflow over the 
northern half of the bay and stronger, more concentrated outflow over the southern half 
of the bay. At depth, geostrophic and ADCF data were consistent on the southern side 
of the canyon and show inflow on the order of 10 cm s -1 . Because of the possible prob¬ 
lems noted in the acquisition of ADCP data along the northern wall of the canyon, the 
geostrophic mean is believed to be more reliable over this area. For this reason, and 
because the geostrophic and ADCP mean data were similar in the other areas of the bay, 
the geostrophic mean field was used to force the numerical ocean model of the Monterey 
Bay. 
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G. NUMERICAL OCEAN MODELING OF THE MONTEREY BAY 
CIRCULATION 


1. Model Description 
a. Dynamics 

The model used in this study is a two layer, semi-implicit, primitive equation 
model which uses a 63 by 37 element C-grid for finite differencing. The numerical 
scheme utilized in this model was initially derived by Hurlburt (1974) to study ocean 
mesoscale circulation features (Bruner, 1988) and later adapted for use in studying cir¬ 
culation processes within the Monterey Bay by Prof. David C. Smith IV of the Naval 
Postgraduate School, Monterey, CA. Linear test cases with analytic solutions have been 
run in tandem by Smith and Reid (1982) to validate the model (Bruner, 1988). The 
equations used in the model are the vertically integrated momentum equation: 


ci 


+ (V- y,+ = + Vi. (12) 


and a continuity equation: 


clti 

cl 


+ V • V, 


0, 


(13) 


for layer (i= 1 upper and i = 2 lower) thickness h„ transports V, and velocities v,. The 
relation: 


P\=g{h\+h 2 +d) (14) 

was used to calculate pressure in the upper layer, and: 

P 2 — P\~ g'h\ (15) 

was used for lower layer pressure calculations where the reduced gravity ( g') is calculated 
using: 


g =fr 


(P 2 -P 1 ) 


Pi 


(16) 


The Coriolis parameter if) is equal to 20 sin </> where Q is the angular velocity of the 
earth in rad s l and <b is the latitude. The Kronccker Delta (()„) is equal to one when i = j 
and zero when /¥= j. Subgrid scale dissipation is modeled using biharmonic (V 4 ) friction 
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with the diffusivity constant ( B „) equal to - 1 x 10 7 m 4 s _1 . Wind stress between the air 
and water (t w ) is included as a forcing parameter along with bottom friction which em¬ 
ploys a linear bottom drag coefficient (r) equal to 5 x 10~ 4 ms _l (Winant and Beardsley, 
1979) in the relation -rv 2 (Shaw and Csanady, 1983). The reduced gravity ( g ') value used 
was 2 x 10‘ 2 mrl 

b. Boundary Conditions and Domain 

The north, east, and south boundaries are closed (no slip). The west 
boundary is open (following Camerlengo and O'Brien, 1980) except where forcing (in¬ 
flow and outflow) is specified. Forcing along the western boundary can be applied sep¬ 
arately to the upper (w'hich encompasses the free surface to 50 m) and lower (50 m to 
the bottom) layers, allowing both barotropic and baroclinic modes. 

The 63 by 37 grid is spaced at 500 m intervals resulting in a model domain 
encompassing a 32 by 18 km rectangle. The domain is orientated roughly north- 
northwest to south-southeast (see Figure 33 on page 62) to encompass as much area 
of the Monterey Bay as possible. Realistic topography (Figure 33 on page 62) was in¬ 
corporated by applying an interpolated field of gridded bathymetry. Because of numer¬ 
ical considerations, the layer interface cannot intersect either the free surface or bottom 
topography. For this reason, bottom topography shallower than 100 m was not allowed 
so unrestricted vertical displacement of the interface could occur (Bruner, 1988). 

2. Initialization Methods 

a. Methods of Model Forcing 

The first problem encountered was whether to use the mean flow obtained 
through geostrophy or through ADCP measurements. Since the sense and magnitude 
of the two means w r ere similar and, because of the doubtful reliability of the ADCP data 
at depth near the northern wall of the canyon (discussed at length in the ADCP section), 
it was decided to use the mean obtained through geostrophy. Aside from that one re¬ 
gion, they were not significantly different. Since the model is a two layer system, the 
mean geostrophic flow was averaged over the top 50 m to obtain v, values, then averaged 
between 50 m and the bottom depth between each station (Figure 34 on page 63) to 
obtain values for v 2 . This averaging was performed in order to conserve the mean mass 
transport of the mean geostrophic flow field. This procedure left the surface velocities 
relatively unscathed, but resulted in lowered velocities at depth since most of the mean 
deep flow lay above 250 m. The next obstacle was how to force the model with known 
data at inconsistent intervals and not in coincidence with the western boundary of the 
model (see Figure 33 on page 62). Since both Line 1 and the western boundary of the 
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Figure 33. Numerical model bathymetry (top) and spatial domain (bottom): Top 
figure shows the bathymetry used in the Monterey Bay numerical 
model, bottom figure shows the domain (box) and its orientation within 
the bay. Line 1 (stations 1-10; and Line 2 (stations 10-17) are also de¬ 
picted in relation to the model domain. Arrows roughly approximate 
the method used in projecting the mean gcostrophic flow onto the 
model's western (open) boundary. 
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Figure 34. Cross sections of model inputs: Top (left) figure shows the result of 
vertically averaging the mean geostropbic flow along Line 1. Top 
(right) figure depicts the projection of these averaged flows onto the 
western boundary of the model. The bottom figure represents the mass 
balanced flow input into the model during subsequent experiments. 
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model are orientated roughly perpendicular to bottom slope, the simplest solution was 
to project the mean (low at known points onto the western boundary of the model grid 
(see Figure 33 on page 62) along isobaths. This of course is a realistic assumption if 
vorticity is conserved over the mouth of the bay and if the predominant flow was per¬ 
pendicular to Line 1. Once the velocities were identified on the model boundary, the 
data were run through a linear interpolation scheme to define velocity values at each grid 
point along the boundary (see Figure 34 on page 63 and Appendix C). Runs were then 
made using this farcing scheme, both with and without the mean wind and bottom fric¬ 
tion (Table 4 on page 65). These model runs ran approximately 8 days before the 
interface intersected the bottom. When this occurred the runs terminated because model 
assumptions were violated. This was not surprising since the original mean geostrophic 
flow was not exactly mass balanced (positive net transport, or flow into the bay), thus 
surplus water was continually being pumped into the bay at each time step with no 
outlet available. In an attempt to stabilize the model, the input mass field (m : s~ l ) ori¬ 
ginally derived from the geostrophic mean flow was adjusted so that net transport in 
both layers added up to zero. This was accomplished by applying linear correction co¬ 
efficients to both layers (Appendix C). This adjustment of mass balance is considered 
justifiable and realistic for the following reasons: 

• The "true" mean flow would be mass balanced. 

• The amount of adjustment applied to the boundary forcing is probably less than 
the geostrophic resolution. 

• As a consequence of the model's two layered system, averaging of the mean flow 
was only conducted down to 50 m for the upper layer. Inspection of the mean 
geostrophic flow field (Figure 24 on page 42) shows that, while it definitely appears 
to be an uncoupled, two layer system, the top layer most likely extends deeper than 
50 m. Therefore, much of the flow in the lower layer was artificially introduced 
during the partitioning process. 

Runs with the balanced flow field appeared more stable numerically, and ran out to an 
average of 19 days. 

To see what the model response would be to wind forcing alone, a run was 
made simulating the actual winds over the bay during a strong frontal passage (similar 
to what was described in the Meteorological Synoptic Situation section). Input winds 
were first "ramped" (i.e. gradually increased each time step to avoid numerical instabil¬ 
ity) from 0 to 2 m r 1 from the south to southeast. Then, from da\ 3 to day 7, the winds 
were ramped from southerly to southwesterly, and the winds were ramped up to 8.5 m 
s~ l to simulate pre-frontal winds. From day 7 to day 8, winds were kept strong, but di- 
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rection was shifted each time step from southwest to northwest to represent frontal 
passage. Winds were gradually decreased from 8.0 mr 1 to 1.0 m s~\ while direction 
was backed to westerly, over the next two days to day 10. Past day 10 the winds were 
turned off so the "spun down" state could be observed. Both layers were initially at rest 
and the entire western boundary was left open (allowing inflow and outflow anywhere). 

All model runs were made with bottom topography included (for 
topographic effects see Bruner, 1988). Open boundary conditions were specified wher¬ 
ever mean forcing was determined to be zero along the boundary to allow inflow or 
outflow at that point. As previously mentioned, the entire boundary was open for the 
wind forcing experiments. The mean MBA wind speed and direction was used in those 
mean forcing runs which applied a constant wind stress throughout. A summary of the 
five runs conducted is provided in Table 4 (outputs appear in Figure 35 on page 66 to 
Figure 39 on page 70). 


Table 4. SUMMARY OF MONTEREY BAY NUMERICAL MODEL SIMU¬ 
LATIONS 


Run No. 

Boundary' Forcing 

Bottom Friction 

Wind Stress 

1 

Geostrophic Mean 

No 

No 

2 

Geostrophic Mean 

Yes 

260T at 2.3 m r 1 

3 

Mass Balanced 
Geostrophic Mean 

No 

No 

4 

Mass Balanced 
Geostrophic Mean 

Yes 

mam 

5 

None 

No 

Variable (see text) 


3. Results of Experiments 

Results of run number one (Figure 35 on page 66) shows anticyclonic rotation 
over the southern half of the bay with a closed cyclonic eddy situated in the northern 
bight region in the surface layer. At depth, velocities are weakly cyclonic (less than 2.5 
cm s _1 ) over the mouth of the canyon with weak anticyclonic curvature in the northern 
and southern bights. The addition of mean winds and bottom friction (run two) did not 
seem to have a marked effect on the overall flow (Figure 36 on page 67). Flow patterns 
w'ere similar in both cases, with the exception that, as expected, velocities w'ere decreased 
in run two by the effect of bottom friction. To gain a quantitative comparison with the 
mean flow' calculated across Line 2, geostrophic velocities w’ere averaged down to 50 m 
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Figure 35. Monterey Bay model output from run No. 1: North is to the left in 
each figure (i.e. figures are orientated as if you are looking into the 
bay). Arrows are velocity vectors, solid lines are isotachs (cm .r 1 ). Top 
figure depicts upper level flow and bottom figure shows deep flow, both 
at day 7 of the run. Mean gcostrophic flow was used to force the model 
along the western (bottom of plots) houndary. The effects of wind and 
bottom friction are neglected. 
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Figure 36. Monterey Bay model output from run No. 2: North is to the left in 
each figure (i.c. figures are orientated as if you are looking into the 
bay). Arrows arc velocity vectors, solid lines are isotachs (cm s -1 ). Top 
figure depicts upper level flow and bottom figure shows deep flow, both 
at day 7.25 of the run. Mean geostrophic flow was used to force the 
model along the western (bottom of plots) boundary. The effects of 
wind and bottom friction arc included. 
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Figure 37. Monterey Bay model output from run No. 3: North is to the left in 
each figure (i.c. figures are orientated as if you are looking into the 
bay). Arrows are velocity vectors, solid lines are isotachs (cm .v 1 ). Top 
figure depicts upper level flow and bottom figure shows deep flow, both 
at day 15.8 of the run. Mass balanced mean gcostrophic flow was used 
to force the mode! along the western (bottom of plots) boundary. The 
effects of wind and bottom friction arc neglected. 
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Figure 38. Monterey Bay model output from run Nr. 4: North is to the left in 
each figure (i.c. figures arc orientated as if you are looking into the 
bav). Arrows arc velocity vectors, solid lines are isotachs (cm •sr 1 )* * op 
figure depicts upper level flow and bottom figure shows deep flow, both 
at day 15.8 of the run. Mass balanced mean gcostrophic flow was used 
to force the model along the western (bottom of plots) boundary. 1 he 
effects of wind and bottom fiiclion arc included. 
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Figure 39. Monterey Bay model output from run No. 5: North is to the left in 
each figure (i.e. figures aic orientated as if you are looking into the 
bay). Arrows are velocity vectors, solid lines are isotachs (cm .r 1 ). Top 
figure depicts upper level flow and bottom figure shows deep flow, both 
at day 13.6 of the run. Wind forcing alone was used to force flow in 
the model. Both layers were initially at rest. The entire western 
boundary' (bottom of plots) was open to allow flow to develop at will. 
The effect of bottom friction was neglected. 










and from 50 m to the bottom and compared with model outputs at these grid points 
(Table 5 on page 71, see also Figure 33 on page 62 for station positions). The mean 
Line 2 data were not used to force the model in any way. The mean geostrophic flow 
between stations 17 and 16 does show light outflow in the extreme northern part of the 
bay, although not to the extent the model shows. Outflow is not indicated in the mean 
Line 1 geostrophic velocities implying that some sort of weak, closed circulation may, 
in fact, exist in the northern bight area. Model output across Line 2 also shows strong 
inflow in the lower layer over the northern part of the bay and outflow in the southern 
part. Line 2 mean velocities show no such flow in these areas. In fact, the mean Line 2 
geostrophic flow shows very little flow in the lower layer at all, with 1.94 cm s~ l inflow 
between stations 15 and 14, and -0.99 cm s _1 outflow between stations 14 and 13. In 
general, the overall flow patterns predicted by the bay model are much stronger, in both 
upper and lower layers, than those found in the observed mean flow field. Circulation 
in the northern bight which was generated in runs 1 and 2 can be explained in terms of 
compression and expansion of vortex tubes. As water is pumped into the northern part 
of the bay (too much water) in the upper layer, it travels into the northern bight region 
where it forces the interface down producing vortex stretching in the upper layer and 
compaction of tubes in the lower layer resulting in cyclonic circulation in layer one, and 
anticyclonic circulation in layer two. 


Table 5. COMPARISON OF MODEL OUTPUT AND LINE 2 GEOSTROPHIC 


MEAN: Negative velocities indicate flow out of the bay, positive values 
denote flow into the bay. Velocities are in cm s~'. _ _ 


Station 

Observed Mean 

Run 2 

Run 4 

Run 5 

v 'i 



V*2 

0 


v. 

Vj 

17-16 

-0.7 

0 

-5 

9 

1 

0 

4 

10 

16-15 

-0.3 

0 

2 

3 

3 

0 

3 

7 

15-14 

4 

2 

10 

0 

4 

0 

1 

-.5 

14-13 

4 

-1 

6 

0 

5 

-.5 

-.6 

-.5 

13-12 

4 

0 

5 

1 

3 

1 

-.5 

-2 

12-11 

2 

0 

4 

1 

1 

-.5 

-1 

12 ! 

11-10 

-9 

0 

-1 

0 

-3 

-.5 

-1 

2 


The eflect of balancing the mass transport along the western boundary 
produced major changes in the flow pattern (Figure 37 on page 68). The circulation in 
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the northern bight produced by run one and run two is no longer evident in runs three 
and four (Figure 37 on page 68 and Figure 38 on page 69). The surface layer consists 
of a broad anticyclone encompassing practically the entire bay, and centered over the 
canyon axis. In run three, flow in the lower layer generally follows topography with 
weak cyclonic circulation over the canyon mouth, somewhat stronger anticylonic rota¬ 
tion over the northern shelf, and anticyclonic circulation contained in the southern bight 
region. Flow appears to diverge from the cyclone in the canyon's mouth, travel deeper 
into the bay along the canyon axis (east), then enters the southern bight via the canyon 
head where it ties into the previously mentioned lower layer anticylone located in this 
area. The addition of bottom friction effectively keeps the flow trapped in the canyon 
(Figure 38 on page 69). Magnitude and direction of flow in both layers from run four 
compares favorably with the observed mean geostrophic flow (Table 5 on page 71) with 
the exception of not predicting weak outflow in layer one over the extreme northern part 
of the bay. The model also appears to under predict the outflow in the southern part 
of the bay in layer one (-3 vice -9 cm s~' which was observed). A Silicon Graphics 
"movie" of this run, made at 1000 s intervals, showed that the model was highly stable 
using this balanced forcing. The predicted high dynamic height center remained quasi¬ 
stationary over the canyon mouth throughout the 19 day animation, agreeing well with 
the inferred flow from the actual mean dynamic height analyses (Figure 28 on page 
49). Again, the model realistically predicts the observed Line 2 mean flow, especially in 
the lower layer, which supports the hypothesis that much of the cyclonic circulation 
within the canyon at depth does not penetrate into the canyon to any appreciable extent, 
but rather recirculates near the mouth of the canyon. 

The final model simulation (Figure 39 on page 70), which used variable wind 
stress as the sole forcing parameter, showed that strong atmospheric frontal passages, 
such as the type previously described, can set up circulation patterns which are similar 
in nature, at least in the upper layer, to those observed in the mean geostrophic flow 
(Figure 39 on page 70). The upper layer outflow aiong the southern portion of the open 
boundary is forced to the southwest, similar to what was depicted in the 04 May satellite 
imagery. Reproducing flow similar to that observed in the lower layer mean geostrophic 
flow using this method was not as successful. The model output shows anticyclonic 
turning in the northern shelf area and cyclonic turning in the lower layer over the 
southern shelf break. These flows then converge near the canyon head and flow out to 
sea along the axis. One major shortcoming of using this forcing method is the fact that 
theie is no model topography above sea level to produce lee effects such as in the 


72 


northern and southern bight regions which, in reality, are fairly protected from 
northwesterly winds (Breaker and Broenkow, 1989). Animation of this simulation 
showed height anomalies to vary widely in location and intensity. This run never quite 
reached a stable state, even after day 10, when wind stress was turned off implying wind 
alone is not the only factor. 





IV. DISCUSSION 


A. COASTAL PROCESSES 

Strong northwesterly winds, produced by post-frontal high pressure ridging across 
California, prevailed over the central California coast from the last week of April 
through the first week of May 1988. The oceanic response to this wind stress was 
coastal upwelling, produced by Ekman transport, and the appearance of a strong 
equatorward geostrophic jet along the central coast. NOAA satellite imagery from 04 
May 1988 clearly shows this upwelling and the presence of equatorward flow, entering 
the Monterey Bay from the north, and advecting warm Bay waters south past Pt. Sur, 
37 km to the south of the Monterey Bay. This pattern is similar to observations made 
by Breaker and Broenkow (1989), who noted that, during upwelling events, cold, up- 
welled water was advected across the bay, producing a strong thermal front near Pt. 
Pinos. A geostrophic velocity cross section, produced from hydrographic data collected 
along a line perpendicular to Pt. Sur on 04 May 1988 depicts a strong 25 cm s~ l 
equatorward flow extending from the coast to 50 km seaward and limited to pressures 
of less than 100 dbar, typical of coastal upwelling geostrophic jets (Figure 40 on page 
75). These descriptions of the winds, upwelling, and coastal geostrophic flow are all 
consistent with the findings and observations of Halliwell and Allen (1987), Buyer 
(1983), and Hickey (1979). Of final note, the CUC was not evident on the 04 May 
geostrophic cross section at Pt. Sur, consistent with Chclton's (1984) observations that 
this poleward, deep current was weak and equatorward during April and May at 
Monterey and Pt. Sur. 

The 08-11 May 1988 Monterey Bay hydrographic survey commenced at the cessa¬ 
tion of strong northwesterly winds. Winds were light and seas were calm throughout the 
period. As a result of this period of calm winds, satellite imagery taken on 08 and 11 
May 1988 showed waters within the bay to have w'armed significantly, either by 
advection of warmer oceanic waters into tne bay (Breaker and Broenkow, 1989), or 
simply through solar radiative processes. While no quantitative evidence exists to sup¬ 
port the continued presence of the equatorward coastal current during this period, weak 
upwelling was still visible to the north and south of the Monterey Bay. 
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Figure 40. Pt. Sur geostropliic velocity cross-section: Geostrophic Velocities were 
derived from C I'D casts collected on 04 and 05 May 1988. The section 


runs roughly east-west, perpendicular to the coast, from Pt. Sur 
seaward. Pt. Sur is located roughly 20 n mi (37 km) south of the 
Monterey Peninsula. Positive values denote poleward flow, negative 
values denote cquatorward How. 


B. INTERNAL WAVES 

The mean flow signal, especial!} over the shelf break and the Monterey Submarine 
Cans on, was highly aliased by the internai wave field throughout the data acquisition 
period (Figure 19 on page 34). 'lire process of averaging the mass field and ADCP ve¬ 
locities appears to have effectively removed the majority of tins 'noise.' Examination 
of the ADCP passes along Line 1 showed that surface tidal heights and the flow within 
the canyon were not in phase. The internal wave amplitudes of 30 m, observed in Line 
1, aic smaller than the 50-120 m amplitudes observed by Shea and Biocnkow (1982) near 






the canyon head. This difference can be explained by their hypothesis that internal wave 
energy is focused near the head as they shoal, thus producing the somewhat higher am¬ 
plitudes they observed. This is also consistent with Broenkow and McKain's (1972) 
observations that internal tidal waves generally decreased seaward. Why observed 
internal wave amplitudes observed in Line 2 data, which is inshore of Line 1, were 
smaller than those observed in Line 1 is unknown. 

C. CIRCULATION 

The mean flow during the data acquisition period is an uncoupled two layer system, 
consisting of a warm core anticyclone at the surface over a cold core, cyclone which is 
"trapped" within the confines of the Monterey Submarine Canyon under the 
thermocline. Inflow at the surface is broad, extending from the northern part of the bay 
south to the center of the canyon, while outflow is relatively stronger but narrower and 
is confined to the area between the Monterey Peninsula and the center of the canyon 
(station 7). Mean flow at depths of 100 to 200 m consists of two cores with outflow 
along the northern wall of the canyon and inflow along the southern half of the canyon 
axis. Both ADCP and geostrophic means (Figure 30 on page 53 and Figure 24 on page 
42), along with dynamic topographies (Figure 28 on page 49), support this description. 
The only major discrepancy noted between the ADCP and the geostrophic means was 
in the outflow region along the north canyon wall. 

Evaluation of the mean flow pattern in Line 2 also shows surface inflow over the 
center of the bay with outflow along the southern portion. Surface velocities are gen¬ 
erally weaker in Line 2 and no evidence of cyclonic circulation at depth is evident, just 
a weak inflow near the center of the canyon. This leads to the conclusion that the sur¬ 
face flow is turning over Line 2, thus accounting for the small east-west velocity com¬ 
ponents, and that the cyclonic flow noted in the canyon itself is limited to the deeper 
portions of the canyon, near the mouth. 

As seen in the cross-sections of mean temperature and density, as well as in the 
mean 10'C isotherm plots, the isotherms and isopycnals had a much more pronounced 
slope from the canyon axis to the Monterey Peninsula explaining the higher geostrophic 
and ADCP velocities found there. Map contours of mean temperature and dynamic 
heights support these flow patterns at the surface and at depth. The tighter dynamic 
height gradient over the southern portion of the bay is also indicative of stronger 
geostrophic flow in this area. Satellite imagery, especially the 04 May 1988 pass, sup 
ports anticyclonic surface flow. The surface flow pattern agrees with the findings of 
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Bigelow and Leslie (1930), Skogsberg (1936), Lammcrs (1971), and with Pirie and 
Stellar's (1974) inferred southward flow from March through July. The noted reversal 
in flow confirms Skogsberg's (1936) suspicions that flows within the bay at different 
depths may be in opposition to one another. Lasley (1977), Smethie (1973), Broenkow 
and Smethie (1978), Reise (1973), Moomy (1973), and Breaker and Broenkow (1989) 
determined flow to be predominantly cyclonic within the bay but most acknowledged 
that periods of flow reversals do occur within the bay. During the period of this study, 
one of these "reversals" was obviously taking place. 

As previously mentioned the CUC was not in evidence, at least at Pt. Sur, during 
the period so momentum transfer from this current does not seem a viable flow gener¬ 
ation mechanism (Breaker and Broenkow, 1989) during the May 1988 cruise. The best 
explanation for the observed flow patterns stems from the results of Klinck's (1989) 
three-layer geostrophic adjustment model of flow over a canyon. Despite the obvious 
differences between modeled shelf and canyon geometry', and the actual bathymetry of 
the Monterey Canyon and vicinity, there are similarities between his results and the flow 
patterns observed in the Monterey Bay. His model results showed that, for a narrow 
canyon (shown to be the case here in Chapter III), anticyclonic surface flow over 
cyclonic, canyon trapped flow would be produced by the geostrophic adjustment process 
that takes place when a surface current encounters a submarine canyon. Similarities 
between the 50 km wide current used in forcing his model, and the 50 km wide, 
equatorvvard surface current depicted in the Pt. Sur geostrophic velocity cross-section 
(Figure 40 on page 75) is also noted. Klinck's results also indicated that once a canyon 
width narrows to less than half of the internal Rossby radius, it has little perturbing ef¬ 
fect on the cross canyon flow. The calculated internal Rossby radius at station 13 was 
11.4 km and the canyon width in this area is approximately 6 to 7 km. This may, in part, 
explain the absence of any significant deep flow at depth in Line 2, where the Monterey 
Submarine Canyon begins to narrow rapidly. 

Whether or not flow reversals within the Monterey Bay are correlated to strong 
upwelling events cannot be determined from the results of this study alone. The findings 
of Broenkow and Smethie (1978) implied predominantly northward flow within the bay 
even during the spring and summer upwelling period. This may suggest that the observed 
mean flow pattern may only develop during periods of strong equatorward, coastal 
flows. The minimum velocity required to drive such flows, the amount of time this 
coastal current requires to "spin-up" the observed flow, or the length of time this flow' 
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pattern remains after cessation of this current, remain unknown quantities and are left 
as topics of future research. 

D. MEAN FLOW MODELING 

The general flow pattern described above was also reflected in numerical model 
simulations, particularly when forced by a balanced flow (runs 3 and 4). The differences 
between the unbalanced (Figure 36 on page 67) and balanced (Figure 38 on page 69) 
runs can be seen by analyzing model outputs of height anomalies (Figure 41 on page 
80). Run 2 height anomalies showed surface heights of approximately 7 cm over the 
mouth of the canyon with a relative minimum of 5.6 cm over the northern bight resulting 
in anticyclonic and cyclonic circulations respectively in the upper layer. The interface 
from run 2 showed a depression of 28 m over the northern bight but 38 m over the 
canyon mouth creating relative highs and lows respectively. Analysis of Run 4 height 
anomalies only showed one high in the upper layer and one low in the lower layer, both 
vertically stacked and situated roughly over the mouth of the Monterey Submarine 
Canyon. These balanced runs seem more realistic than those produced by the unbal¬ 
anced mean forcing (runs 1 and 2) and depict broad anticyclonic flow in the upper layer 
which effectively fills the bay. These runs also appear to agree better with observed mean 
flows derived from Line 2, and also with geostrophic flows implied from dynamic to¬ 
pographies (Figure 28 on page 49). The strong southwesterly outflow near Pt. Pinos, 
noted in the observed mean flow and in satellite imagery, appears to have been handled 
well by the model. When the effects of bottom friction were added, the majority of the 
deep flow was confined to the canyon topography, also supported by the observed weak 
flow at depth in the Line 2 mean. The effects of friction are similar to those noted by 
Garcia (1971) who observed that in deep areas of the bay, friction played a small role 
but in shallower areas, friction significantly altered the flow. Some evidence exists in the 
mean Line 2 data to support outflow (or return flow) along the northern coast of the 
bay where the inflow diverges however not to the extent of that predicted by the model 
from the unbalanced mean forcing (runs 1 and 2). 

The major limitations of these model results include: 

• The assumptions used in projecting the mean flow field onto the model's open 
boundary. The mean flow inferred by Line 1 data (west of the open boundary) 
may not, in reality, have extended all the way into the model's western boundary, 
as was assumed when forcing was applied. However, in view of satellite imagery 
and dynamic height analyses, it is felt that model flow was fairly representative of 
the actual flow encountered over the period of data acquisition. 
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• The upper layer actually extends to the bottom over the shelf regions of the 
Monterey Bay however, due to model assumptions, these areas must be treated as 
two layers to avoid surfacing or grounding of the layer interface. 

• The observed flows indicate that a three-layer model is actually required, i.e. the 
top layer to model flow over the shelf area, a second layer to model flow at inter¬ 
mediate depths of 100 to 500 or 600 m, and a bottom layer to model flow in the 
quiescent regions of the Monterey Submarine Canyon below about 600 m. 

Results of these runs are in general agreement with Bruner's (1988) observations that, 
when lower layer velocities are weak and vertical shear is large, the model effectively 
"decouples" the two layers, resulting in independent flow patterns. Bruner noted that, 
for narrow inflow to the north, and broadened outflow to the south, anticyclonic circu¬ 
lation ensued over the entire bay. In contrast, anticyclonic flow also resulted in the runs 
produced in this study for broad inflow over the northern bay, and narrow outflow over 
the southern bay. 

Finally, model runs which were forced solely by variable wind stress, while obtaining 
similar mean flow patterns in the upper level, did not correspond to observed mean flow 
at depth. The wind forcing did, however, show decoupling of the two layers. Obviously, 
other factors besides wind stress were at work here but these experiments do imply that 
wind forcing may play some part in shaping mean flow patterns, at least at the surface, 
during strong wind events lasting on the order of a week. Another interesting aspect of 
this model result is the fact that approximately 4 days after the wind forcing was turned 
off (Figure 39 on page 70), flow remained relatively strong in both layers implying th 
"signal" of these wind events may be seen in the mean flow for a reasonable amount of 
time after the winds have abated. Modeled current flow' appeared to respond rapidly to 
changes in both wind direction and magnitude, agreeing with Hickey's (1979) observa¬ 
tions that oceanic responses to changes in wind stress are rapid (O(days)). However, as 
noted in Chapter III, model runs using only wind stress forcing did not achieve any ap¬ 
preciable level of stability and displayed highly variable height anomaly patterns. 

From modeling results it can be seen that the combined effects of winds, bottom 
friction, and, of course, boundary forcing all play key roles in producing mean circu¬ 
lation patterns within the Monterey Bay. Strong winds appear to have a significant ef¬ 
fect on surface flow, while light mcsoscale land and sea breezes typically encountered 
under weak high atmospheric pressure gradients have little effect on surface flow. Bot¬ 
tom friction played a key role in suppressing lower level flow out of the Mon^rey Sub¬ 
marine Canyon. It was also shown that the model was highly sensitive to mass balanced 
flow, yielding strikingly different flow patterns from balanced and unbalanced flows. 
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Figure 41. Run 2 (left) and run 4 (right) height anomalies.: Left figuics depict 
Run 2 height anomalies (m) for the upper layer (top) and bottom layer 
(bottom). Right figures show Run 4 height anomalies for upper layer 
(top) and bottom layer (bottom). 


Finally, while the model used in this study could not be forced by an alongshore current, 
this component was probably at woik during the period of this study. Actual flow was 
undoubtably influenced by an alongshore geostrophic jet and, as previously stated, sim¬ 
ilarities to results obtained by Klinck (1989) arc noted. This near coastal jet was prob¬ 
ably the major generation mechanism of the observed mean flow patterns. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Both prior to and during the period studied in this report, upwelling was occurring 
along the entire central California coast. Both satellite imagery and the 04 May 1988 
geostrophic velocity cross section at Pt. Sur supports this conclusion. Within the 
Monterey Bay itself, the observed internal wave field was quite active and effectively 
masked the mean signal throughout the period. Previous studies suggest these internal 
waves are quite common within the bay and not an anomalous situation during this 
study. Averaging the mass field and measured ADCP flows effectively filtered out this 
"noise." Because of this high noise level, any type of hydrographic data acquired over 
a short time period from the bay is of little use. Only long term averaged data can be 
used to infer any type of semi-permanent flow pattern w'ithin the bay. Also, because of 
the large amplitudes of these internal waves, geostrophic velocities derived from single 
section mass field observations will typically yield unrealistically high values, especially 
within the thermoclinc. 

During periods of strong surface inflow in the northern part of the bay, the flow 
within the Monterey Bay becomes a two-layered uncoupled system consisting of 
anticyclonic How at the surface and weak, canyon trapped, cyclonic flow below the 
thermocline, limited mainly to the deep waters over the canyon mouth. The main 
mechanism responsible for this flow pattern during May 1988 was the alongshore, 
geostrophic, equatorward flow turning into the northern half of the bay. As this 
equatorward flow encountered the Monterey Submarine Canyon, the geostrophic ad¬ 
justment process caused upwelling within the canyon which resulted in vortex tube 
compaction (anticyclonic flow) in the upper layers and stretching of vortex tubes at 
depth (cyclonic flow) below the thermocline (Klinck, 1989). Numerical simulations of 
the mean flow pattern suggest that the effects of bottom friction may act to confine 
lower layer flows to the canyon itself. Other numerical simulations suggested that wind 
stress may induce anticyclonic surface flow when pre-frontal and post-frontal winds were 
used to force the model although these experiments were inconclusive due to lack of 
evolution to a steady state. 

If surface flow within the Monterey Bay is predominantly cyclonic as previous au¬ 
thors suggest (Breaker and Brocnkow, 1989), upwelling periods nav coincide with p i- 
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ods of flow reversals that have been observed in previous studies (Breaker and 
Broenkow, 1989). The amount of time this flow reversal lasts after the cessation of up- 
welling and the alongshore geostrophic jet is unknown. The suggested scenario leading 
up to the mean flow pattern encountered in this study involves: 

• A 7 to 10 day period of strong northwesterly winds produced by high pressure 
ridging across California behind a series of atmospheric cold fronts. 

• Wind stress induced Ekman transport resulting in strong upwelling along the cen¬ 
tral California coast. 

• Development of a strong geostrophic, equatorward, coastal jet along the central 
coast. 

• A geostrophic adjustment process, such as that described by Klinck (1989), as this 
coastal jet encounters the Monterey Submarine Canyon, resulting in the flow pat¬ 
terns described above. 

Since winds are upwelling favorable during most of the summer season, anticyclonic 
circulation within the Monterey Bay may be more common than previously believed. 

B. RECOMMENDATIONS FOR FUTURE STUDIES 
1. Currents 

The internal wave field within the Monterey Bay needs to be resolved. Once a 
solid grasp of the frequency and wavenumber spectrum is known, different sampling 
schemes can more effectively cancel out their influence on the mean field. In this way 
mean flow can be resolved with a minimum number of samples, thereby saving expensive 
ship time. To accomplish this, an array of fast sampling thermistor chains should be 
deployed at strategic locations throughout the bay such as on the shelf breaks. An al¬ 
ternative approach to averaging out the internal wave aliasing is simply through statis¬ 
tical averaging. A higher sampling population should result in smaller standard errors 
and a better resolution of the mean flow, although cost constraints could be more of a 
consideration in this case. 

Because of non-linear processes at work within the Monterey Bay, current me¬ 
ters and ADCP velocities are probably of somewhat more value than geostrophy in de¬ 
termining flow within the bay. The better resolution of the mean flow by ADCP 
techniques was apparent in the smaller calculated standard errors of this data set. 
Again, the longer the time series, the better. Bottom moored ADCP data, especially 
along the shelf break, would yield valuable information concerning the cross-canyon 
velocities and upwelling, if any, that is taking place from canyon depths. Seasonal 
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spectral analyses of the MBA time series data should continue to reveal much about the 
shelf processes that are occurring within the Monterey Bay. 

Finally, as with any scientific hypothesis, results should be reproducible. 
Therefore, a similar study, conducted during or shortly after a comparable upwelling 
event, should be undertaken to see if flow patterns similar to those found in this study 
exist. In this way, flow "reversals" within the Monterey Bay could perhaps be correlated 
with strong upwelling events. Likewise, a study conducted during a strong poleward 
flow event would also be enlightening. Current meter moorings placed on the shelf re¬ 
gions to the north and the south of the Monterey Bay during the period of these studies 
could verify correlations, if any, between these alongshore flows and circulation within 
the bay. 

2. Numerical Studies 

While many insights into the behavior of flow within the Monterey Bay were 
gained from the numerical model used in this study, a few added features would greatly 
enhance its meaningfulness, namely: 

• Future studies which use actual flow to force this model should be from data along 
a line coincident with the model's open boundary. This would negate the necessary 
assumptions made in this study in projecting the flow field onto the model grid. 

• As mentioned earlier, it would appear that a three-level model would be more ap¬ 
propriate in modeling the Monterey Bay currents. This would diminish the effects 
of vertical averaging and keep flows closer to their "true" velocities. 

• The addition of a realistic coastline and surface topography would enhance the ef¬ 
fectiveness of wind forced flow experiments within the bay. 

• Since local heating may play a role in circulation processes within the bay (Bigelow 
and Leslie, 1930), the addition of thermodynamic equations to model air-sea heat 
exchanges would most likely yield interesting results. 

• Expanding the model domain used in this study westward, to encompass the deeper 
seaward portions of the Monterey Bay Submarine Canvon, and northward and 
southward, to model the shelf regions in these areas, would provide more realistic 
results. The additional capability of forcing the model with northward or south¬ 
ward alongshore currents, similar to Klinck's (1989), would also be of great value. 
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APPENDIX A. MEAN GEOSTROPHIC VELOCITIES 


The following pages contain the mean Gcostrophic velocities for the section across 
the mouth of the Bay (stations 1 to 10). Velocities were cleaved by calculating the mean 
mass field (mean temperatures and salinities) from data collected during the 08-11 May 
19S8 NTS Student Cruise, and then using the thermal wind relation to calculate the ve¬ 
locities. 
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-10.Cl 

-0.620*5E*04 

-0.50;«E*05 


coo. 

0.0*58 

e.ocse 

0 5?0o1 

0 56842 


-10.40 

- 0. 62t>°5£*04 

-0.4403 7£*05 


CCS. 

O.CCSC 

O.CCSC 

C.53481 

0.532*2 


-io.ee 

-0 6345?E*04 

-0.37768E*05 


250 

o.c :<o 

0.034* 

C.i°°o6 

0.4«?42 


-10 65 

-0.63716E*04 

-o.ji«::e.o5 


C 75. 

o o:<i 

0.032° 

0.4650b 

0.4*285 


-10.51 

-0.606°CE*04 

-0.C5051E*05 


CCC. 

0.0334 

0.0CC5 

C.420°4 

0.428°) 


-° 65 

-0.56C6CE*04 

-0.18«8:E.05 


CCS 

o.c::° 

0-0CC8 

o.:°7i° 

0. C°55° 


-°. 03 

-0 5C°01E*04 

-0.1225oE*05 


:so. 

0.0222 

0 0322' 

0.2*458 

0.2b2 7° 


-8.52 

-0.50«3oE*04 

-0.8065$E*0* 


■ ;c _ 

o ecu 

o.cci; 

0.CCC21 

0.CCC55 


-8.38 

-C.50875E*04 

-0.C°71°E*04 


4C0 

C.OCll 

c.cco? 

0.30Gt8 

0.C°08° 


-8.51 

-C.4SC81E•04 

0.31155E*04 


425. 

o o:ci 

0 CC°8 

0 2t°56 

0.2*01° 



-0.310°CE*C4 

0.66437E*04 


450. 

0 02°7 

0.02°I 

0. CC°15 

0.23825 


-2 80 

-C.14C71E*04 

0.«75C8E*04 


«7S. 

C.C2&° 

o oces 

0 20°<l 

O.CO°Cl 


-0. °4 

0.60063E*OC 

0.111801*05 


SCO. 

O.OC81 

O.OC78 

C. 1804° 

C 18073 


1.14 

0.10°C4£*04 

0.111CCE*05 


SC5 

c.ccrc' 

0 026° 

0.15244 

o is:°o 


C 4° 

0.18°61E*04 

0.100C8E*05 


550 

0.02*4 

0.02b3 

0. 1 CSC 7 

0.12*07 


2.81 

o. cscoje*04 

0 8I315E*04 


575. 

o.ccss 

0.CCS 7 

0.0°885 

0.0°°21 


4 56 

0 C4561E‘C4 

0.561ICE‘04 


cCO. 

o.oco 

O.CCSC 

0.07334 

0.0740° 


2 5° 

0.17l°7E*04 

0 .31551E*04 


bCS. 

0.0245 

0.C247 

0.0C8-S 

0.0488° 


C.1C 

0.10C6CE*04 

0.14254E*04 


tSO 

0.023° 

0 CCiC 

G.023°4 

0.02422 


1.22 

0.3°°1CE*03 

0.3°°1CE*03 


6 75. 

O.00C8 

0.0236 

o.oooco 

O.OOOCO 


0.00 

O.OOOOCE*00 

0.00000E*00 

0 





TOTAL 

TRANSPORT * -0. 

70078E*u5 


c 

VELOCITY C/.LCULA; 

’ED ASSUMING 

A LEVEL OF NO 

MOTION O p t» 7 5 - 

CB 





d: 

IST/NCE EETUEE’I 

STATIONS • 

2.40° KM 







c 


"* 0E0CTRPH1C VELOCITY 

BETWEEN STATIONS 6 AND 7 




o ricei 

AV DV*i( 

6) AV DYNI 7) 

SUM DYNI 

6’ 

SUM pvM 7) 

VELOCITYICM/S) 

TRANSPORT(M3/S) 

SUM OF TRANS 

c 

C. 058° 

C-CbCO 

1.2152° 


1.22080 

10.47 

C 13°4°E*C5 

O.747C0E*O5 

:s 

0. CSC 7 

0 0517 

1 . ISbAb 


1 16077 

8 24 

0.10°27E* 05 

0.40;?IE>05 

50. 

0 0450 

0.045: 

1 .10576 


1.10°C7 

6 .<• C 

0.°0085E *04 

0.49844E«05 

75 

0 041° 

0.C4CC 

1 .0607° 


1 0637b 

5.75 

0.80°36E*C4 

O.<0aJ6E*05 

100. 

0.03 °6 

0 C40C 

1.018°: 


1.0215° 

5.16 

0 67187E*04 

0.2C74CE*05 

ICS 

0 0383 

0 03°0 

0 °7°36 


0. °813 7 

2. °0 

0.468b4E*C4 

0.260C3E*05 

150 

0.0372" 

0.037° 

0.°410° 


0 °42C4 

2.42 

0.27117£*04 

0.Cl 22 7E*05 

175 

0.02*5 

0 0371 

0. °C30S 


0 °C44 ° 

1 24 

C 97184E*03 

0.18625E*05 

COO. 

0 0358 

0 O3o3 

0 8*737 


0 8b/4J 

0.03 

-0.5 ?250E*03 

0.17654E*C5 

CCS 

0.0 352' 

o o:s6 

0 8 215 7, 


0 83112 

-0 C5 

-0.17°14E*04 

0.18226E*05 

CSC 

0 0346 

0 0250 

0 7 ° 6 3 7, 


0 ?°5So 

-l 5 7 

-0 28°17E *04 

0.20017E*05 

C75 

0 033° 

0.0341 

0 7617° 


0.7605° 

-2 23 

-0.2?221E*04 

0.22°0°E*C5 

200. 

0.02:3 

0 0 3 34 

0.72786 


0.72646 

-2 70 

- 0.408 7 °E * 04 

0.2664 )£ *05 

3 C S. 

0 0328 

0.0325 

0.69452 


0.69309 

-2.81 

-0 28007E*04 

0.3072°E*0S 

350 

0 occc 

0.0316 

0.66172 


0 66052 

-2.21 

-0.256685*04 

0.34530E*05 

2 75 

0.0317 

0.0310 

0 62°4° 


0.628°0 

-1.15 

-0.68 2l3E•03 

0.270°7E*05 

400 

0.0307 

0 0301 

0 5°733 


0.5°7°5 

0.22 

0.11C30E*04 

0.37780E*05 

425 

0.02°8 

0.02°5 

0.56714 


0.56780 

1 2° 

C.C2238E *04 

0.36657E*C5 

450. 

o o:°i 

o c:°o 

0.5372° 


0 5292b 

1 fib 

0 C°°44E*04 

0.242C3E*05 

4 75 

0.0285 

0 C28C 

0 5C816 


0.50°28 

2 18 

0.36C25E *04 

0.312C°E*05 

500. 

0.C2 78 

C 02 7 7 

0 47°c8 


0.481C8 

2.71 

0.41080£*C4 

0.27705E»05 

525. 

0.026° 

0.0271 

0.451°1 


0.45237 

2 83 

0.3841SE•04 

0.235°7E*05 

550. 

0.026 3 

0 CCoS 

0 42501 


0.42*22 

2 25 

0.3101CE*04 

0.19?56E*05 

5 75 

0.C257 

0 02 bO 

0.2°875 


0.2°??0 

1 84 

0.C25b 7E *04 

0.16654E*05 

6C0. 

0 0252 

0.0254 

0 27204 


0.2727J 

1 34 

0 !?696£*04 

0.14C°CE*C5 

6C5 

0 024? 

0.0248 

v.24784 


0.24828 

1.04 

0.12887£*04 

0.1C52°E«05 

650. 

o.c:-.: 

0.0241 

0 22316 


0 3235° 

0.83 

0.1417 °E * 04 

0.11140E *05 

6 75 

0 CC3t> 

0 0235 

0 2° Q °4 


0.2°°50 

t .C8 

C 1 682 IE *04 

0.97CC4E»C4 

7C0 

C 0231 

o c::i 

0 27538 


C “75°° 

1.1° 

0 17 * 02E * 04 

0.802°2E*04 

725 

0 CCC? 

0 022 7 

0.25227 


0 25288 

1 .1° 

o.i*;°i£*04 

0.627°1E*04 


S5 








7S0. 

0.0223 

0.0224 

0.22°58 

0.2201A 

1.0b 

0. 14688E*C4 

0.A6000E * 04 

r;s. 

0.021° 

0.0221 

0.20728 

0.20776 

0. °2 

0.UC74E*04 

0.211 ICE * OA 

800. 

0.0217 

0.021° 

0.18525 

0.18564 

0.57 

0.651°5E *02 

0.20038E*04 

b:s. 

0.0C1S 

0.0215 

0.16262 

0.1627° 

0.31 

T.3°A8AE*02 

0,125 1 °E*0A 

8S0. 

0.0211 

0.0212 

0.1*212 

0.14224 

0.22 

0 lb362E*02 

0.95702E*03 

875. 

0.020° 

0.020° 

0.120°° 

0.1209<- 

0,00 

0.2A526£*02 

0.792A1E*03 

“CO. 

0.0206 

0.0205 

0.10011 

0.10013 

0.03 

0 . °8°1OE *02 

0.76888E*03 

°:s. 

0.0203 

0.0202 

0.07°54 

0.0795° 

0,1° 

0.20257E*03 

0.669°7E*03 

°50. 

0.0200 

0.0200 

0.05928 

0.C5937 

0.17 

0.22347E*03 

0,46760E*0J 

0 75. 

0.01*8 

0.01 °8 

0.02927 

0.03°33 

0.13 

0 . 1 707 1E*02 

0.2AA15E*03 

1000 

0.01 °5 

0.01 °5 

0.01*49 

0.01°54 

0.10 

0.7341°E * 02 

0.7341°E*02 

ic:s. 

0.0068 

0.01 °3 

0.00000 

0.00000 

0.00 

o.oooooe*oo 

0.OOOOOE^OO 


0 TOTAL TRANSPORT s 0.74?20E*C5 

0 VELOCITY CALCULATED ASSUMING A LEVEL Of NO MOTION Of IOCS. DB 
DISTANCE BETWEEN STATIONS : 5.932 KM 


0 **••* GEOSTRPMIC VELOCITY BETWEEN STATIONS 7 AND 8 


0 :(DB) 

AV DYN( 7 ) 

AV DYN( 8) 

SUM DYNC 7) 

SUM DYNt 8) 

VELOCITY(CM/S) 

TRANSPORT(M2/S) 

SUM OF TRANS 

0. 

0.0600 

0.058A 

1.0°°81 

1.0°5°8 

-11.12 

-0.8722°E*04 

0.101S7E * 06 

25. 

0.0517 

0.0A°0 

1.03978 

1.02752 

-6.54 

-0.2*44°E*04 

0.!1011E+06 

50. 

O.OA52 

0.0A27 

0.988C8 

0.°8849 

1 .1° 

0.35282E *04 

0.11275E*06 

75. 

0.0A22 

0.0412 

0.°A277 

0. °4482 

5. °5 

0.7320CE*04 

0.10922E*06 

100. 

0.0A02 

0.0400 

C.90061 

0. °C3*t- 

8.8* 

0.9I4J1E*04 

0.101°0E *06 

125 

o.o:«o 

0.0289 

0.86029 

0.86271 

0 .64 

0.97789E*04 

0.92?61E*05 

15C. 

o.o;?° 

0.0290 

0.82125 

0.82485 

10.15 

0. °887°E*04 

0.82°82E*05 

175. 

0.0271 

0.0272 

0.78250 

0.7868° 

9.86 

0.92898E*04 

0.73094E*OS 

2CC. 

C. 0262 

0.03*7 

0.74642 

0.74956 

9.14 

0.84024E*04 

0.63705E*05 

225. 

0.0256 

0.02S8 

0.71014 

0.712§5 

7.86 

0. ?484°E*04 

0.55302E*05 

CSC. 

0.0250 

0.03A° 

0.67457 

0.67708 

7.28 

0.725I2E* 04 

Q.4781?E*05 

275. 

O.02A1 

0.03A3 

0.63°b0 

0.64215 

7.39 

0.70117E*04 

0.40564E«05 

:cc. 

0.022a 

0.0227 

0.60547 

0.60781 

6.7° 

0.62234E*04 

0.S5553E*05 

225. 

0.0225 

O.C22? 

0.57209 

0.S74U 

6.00 

0.56917E*04 

0.2722°E*05 

250. 

O.C216 

0.0219 

0.53954 

0.54144 

5.52 

0.51228E-04 

0.21537E*05 

S7S. 

0.0310 

0.0211 

0.50791 

0-50°58 

4.85 

0.45097E‘04 

0.16415E*05 

ACC. 

0.0301 

0.020A 

0.476°6 

0.47843 

4.28 

0.3°082E*04 

0.1 1905E*05 

a:$. 

C. C2°5 

0.02°7 

0.44681 

0.44806 

3.62 

0.2302SE*04 

0. 79°6?E*04 

ASO. 

o.o:°o 

0.C2° 1 

0.41727 

0.41822 

3.05 

0.28108E*04 

0.46942E*04 

A75. 

0.C282 

0.C2S5 

0.2892° 

0.38°20 

2.64 

0 21792E*04 

0.18834E*04 

see. 

0 0277 

0.027° 

0.2600° 

0.3t>C?0 

1.77 

0.14765E*04 

-0.29S8SE*05 

525. 

0.0271 

C.C27A 

0.23228 

0.52280 

1.22 

0.°0058E*03 

-0.17?24E*04 

550. 

0.02*5 

0.02*b 

0.20S24 

0.20S4S 

0.61 

0 47044E«03 

-0.26?2«E*04 

575. 

0.02*0 

0.0260 

0.27871 

0.27892 

0.25 

0.35?66E*03 

-0.31424E*04 

600. 

0.025A 

0.0252 

0.25274 

0.25287 

0.28 

0.38705E*03 

-0.3$01QE*04 

625. 

0.02A8 

C.02A3 

0.2272° 

0.22752 

0 41 

0.3246°E*03 

-0.38881E*04 

650. 

0.02A1 

O.02A2 

0.202*0 

0.2C270 

0.27 

-0.8*020E*02 

-0.42228E*04 

6 75 

0.0225 

0.0227 

0.17851 

0.1782b 

-0 AS 

-0.74780E*03 

-0.413676*04 

TOO. 

0.0231 

0.0221 

0.15501 

0.15464 

-1 v7 

-0.10707E*04 

-0.32890E*04 

725. 

0.C227 

0.0226 

0.1219° 

0.13151 

-1.10 

-0.°2449E*03 

-0.231836*04 

750. 

0.022A 

0.0223 

0.10° 15 

0.1098° 

-0 77 

-0.6$159E*03 

-0.139386*04 

775. 

0.0221 

0.0220 

0.08677 

0.09658 

-0.55 

-0.388S0E*03 

-0.742216*03 

800. 

0.021° 

0.0218 

0.064*5 

0.06457 

-0.24 

-0.19064£*03 

-0.35391E*03 

825. 

0.0215 

0.0215 

0.04280 

0.04275 

-0 15 

-0.1I804E*03 

-0.16327E*03 

850. 

0.0212 

0.0212 

0.02125 

0.02122 

-0 0° 

-0.45235E*02 

-0.45235E*02 

875. 

0.020° 

o 

o 

l> 

o 

CO 

0.00000 

0.00000 

0 00 

0.OOOOOE*OQ 

0.00000£*00 


TOTAL TRANSPORT * 0.10127E»06 


0 VELOCITY CALCULATED ASSUMING A LEVEL OF NO MOTION OF 87S. DB 
DISTANCE BETWEEN ST AT 10..S : 3.°S 3 KM 


0 



GEOSTRPMIC 

VELOCITY 

BETWEEN : 

STATIONS 8 AND 9 




0 

2(DP) 

AV 0YN( 8) 

AV OVN< 

9) 

SUM DVN( 

8) SUM DYN( 9) 

VELOCITYfCM/S) 

TRANSPORT(M3/S) 

SUM of tran: 


0 . 

0 0584 

0.0S69 


0.81-MS 

0 .61434 

- 12.04 

-0 585516*04 

0.9S366E*04 


25. 

0.04°0 

0.0482 


0.75970 

0.75743 

-5.45 

-0.24465E*04 

0.15192E*05 


50. 

0.0427 

0 0421 


0.70°66 

0.70923 

-1 86 

-0.3°201E*C3 

0.1HJ8E*05 


7$. 

0.C412 

0.04C9 


0-665°° 

0 66*15 

0 63 

0.10248E»04 

o.180S0E*05 


80 








100. 

O.CAOO 

0.03*8 

0.^2683 

0.62528 

2.28 

0.18674E*04 

0.17005E-C5 

ICS. 

0.028? 

0.C3®0 

0.58498 

0.535*2 

2.20 

0.1 ?83?E*04 

0. 15128E*05 

150. 

0.0280 

0 028C 

0.54t0C 

0.546*5 

2.72 

0.14*J8E*04 

0.13154E*0S 

175. 

0.0372 

0.027C 

0.5030* 

O.5C065 

1.65 

0.12794E*04 

0.116®1E*05 

COO. 

0.0267 

0.0263 

0.47074 

0.67124 

2.17 

0.20578E*04 

0.10411E*05 

CCS. 

0.0358 

0.0256 

0.42402 

0.43495 

3.*7 

0.28«41E*04 

0.83534E*04 

C50. 

0.024? 

O.O2S0 

0.2*825 

0.39*34 

6 67 

0.20612E-04 

0.545?3E*04 

C75. 

0.024J 

0.026' 

0.26222 

0.26426 

6.47 

0.2®«S*E‘04 

0.23<7?E*04 

200. 

0.0227 

0.0228 

0.228*8 

0.23002 

4.47 

0.278®2E*04 

-0.597*2E 4 03 

2C5 - 

0.02C7 

0.032® 

0.2*532 

0.29*23 

2.86 

0.23172E*04 

-0.3386*E*04 

:so. 

0.0319 

0.02C1 

0.C6C61 

0.26232 

3. OC 

0.172®8£»04 

-0.57041E*04 

275. 

0.0211 

0.0314 

0.23075 

0.22125 

2.13 

0.1106?E*06 

-0.744J?E*04 

4C0. 

0.0206 

0.0207 

0.1**60 

0.1**68 

1.17 

0.22530E*C3 

-0.85508E*04 

ACS . 

O.OC97 

0.0202 

0.1**24 

0.16*1* 

-0 20 

-0.72630E‘03 

-0.88761E*04 

450. 

0.02*1 

0.02*4 

0.12*40 

0.12*02 

-2.00 

-0.17224E*04 

-0.813«8E*04 

6 75. 

O.OC85 

0.0C8& 

0.11037 

0.10*64 

-3.15 

-0.21831E*04 

-0.64164E*04 

500. 

O.OC7? 

0.0C78 

0.08187 

0.08108 

-2.37 

-0.21536E*04 

-0.42234E*04 

525 • 

0.0C74 

0.0C70 

0.052*7 

0.05326 

-3.06 

-0.1S52*E*04 

-0.20798E 4 04 

sso. 

0.0C66 

O.OCoC 

O.OCooC 

0.02625 

-1.57 

-0.52 724E *03 

-0.52724E*03 

575. 

O.OCoO 

0.0C54 

C.CCOCO 

0 00000 

0.00 

O.OOOOOE*00 

0.0000OE*O0 





TOTAL 

TRANSPORT s 0. 

?3366E*04 



0 VELOCITY CALCULATED assuming a level cf no motion of STS. DB 
DISTANCE BETWEEN STATIONS • 2.678 KM 


0 GECSTRPMIC VELOCITY BETWEEN STATIONS * AND 10 


0 Z(C?) 

AV DVN{ ?) 

AV CVN< 10) 

SUM DYN< ?) 

SUM DYMC 10) 

VELOCITY(CM/S) 

TRANSPORT(M2/S) 

SUM OF TRANS 

0. 

0.C55® 

0.0535 

0.1481? 

0.12**® 

-18.57 

-0 1*4?*E*05 

-0.313*8E 4 05 

25- 

0.0482' 

0.0647 

0.0*123 

0.08*1® 

-11.11 

-0.95?85E*04 

-0.11900E*05 

50. 

0.0431 

0.04 15 

© 

o 

o 

C.04147 

-3.50 

-0.22011E*04 

-0.22011E*04 

75. 

0.04C8 

0.0141 

0 00000 

0.00000 

0.00 

O.OOCOOE-OO 

0.00000£ 4 00 

0 




TOTAL 

TRANSPORT * -0.2] 

[298E*C5 



0 VELOCITY CALCULATED ASSUMING A LEVEL 0* NO MOTION OF ?S. DB 


87 



The following pages contain the mean Gcostrophic velocities for the inland section, 
towards the head of the Monterey Submarine Canyon (stations 10 through 17). Veloci¬ 
ties were derived by calculating the mean mass field (mean temperatures and salinities) 
from data collected during the OS-11 May 1988 NPS Student Cruise, and then using the 
thermal wind relation to calculate the velocities. 


DISTANCE BETWEEN STATIONS r 5.256 KH 


l >--»• CEOSTROPMIC VELOCITY AND TRANSPORT FOR DATA OF THE CRUISE MAY 1989 

(UNITS IN DYNAMIC METERS) 

o mm.mm OEOSTRPHlC VELOCITY BETWEEN STATIONS 17 AND i6 mmmmm 


0 Z(DP) 

AV DVM( 1/) 

AV D\M( 16) 

SUM DVN< 

17) 

SUM DVNC 16) 

VELOCITVCCM/S 

) TRANSPORT(MI/S) 

SUM OF TRANS 

0 . 

0 OSo 1 

0.0554 

0.05611 


0.05525 

-2 04 

-0.10851E*04 

-0.10851 E*04 

25 . 

0 0CS8 

0.0457 

0.00000 


0.00000 

0.00 

0.0C000E *00 

0.OOOOOE*QO 

C 





TOTAL 

TRANSPORT » -0. 

1085 IE* 04 


0 VE 

LOCI TV CALCULATED ASSUMING 

A LEVEL OF NO 

MOTION OF 25. DB 




01 

STANCE between 

STATIONS s 

4.246 KM 















o r(tp) 

AV DYN( 16) 

AV DYN( 15) 

SUM DYN( 

16) 

SUM DYM( 15) 

VELOCITY(CM/S 

) TRANSPORT(MJ/S ) 

SUM OF TRANS 

0 . 

C.CS54 

0.05*4 

0.10101 


0.10054 

-1.55 

-0.4°84JE«05 

-0.S32°8E*02 

25. 

0.0*57 

0.0460 

0.04566 


0 . 045*3 

0.75 

0.4b513E*03 

0.465136*03 

5C. 

0.0U6 

0.0415 

o.oocoo 


0.00000 

0 00 

0.000006*00 

0.000006*00 

0 





TOTAL 

TRANSPORT t -0. 

332°8E*02 



0 VELOCITY CALCULATED ASSUMING A LEVEL OF NO MOTION OF 50. DB 


DISTANCE BETWEEN STATIONS : 4.9t>3 KM 


c 

Ml 

GEOSTPPHIC 

VELOCITY 

BETWEEN STATK 

5N$ 15 AND 14 

N «* •> * 



0 2(DB) 

AV D>«( 

15) AV DYN< 

14) 

SuM DVM( 15) 

SUM DYN( 14) 

VEL0CITY(CM/3) 

TRANSPORT(M3/S) 

SUM OF tran: 

0 

0 0544 

0.0531 


0. 43°00 

0.43°22 

0.77 

0.24230E *C4 

0.236°5E*05 

25. 

0.0460 

0.0457 


0.38465 

0.58613 

5.22' 

0.4605*E *04 

0.212726*05 

50. 

0 0415 

0.0418 


0 338*7 

0.34041 

6.16 

0 4563*E*04 

0.166666*05 

75. 

0.C3°l 

0.03«7 


0.29717 

0.2°862 

5.12 

0.32365E *C4 

0.121026*05 

100. 

0 0384 

0.038* 


0.25800 

0.258°l 

2 89 

0.20456E*C4 

0.88659E‘04 

125. 

0 C379 

0.0378 


0.21071 

0.22033 

2.17 

0.I°4«4E*04 

0.682036*04 

150. 

0 0373 

0.0373 


0.18183 

0.18258 

2.64 

0.200476*04 

0.487096*04 

175. 

0.03*7 

0.0370 


0.14457 

0.14523 

2.31 

0.154966*04 

0.286626*04 

200. 

0.0364 

0.0566 


0.10784 

0.10827 

1.52 

0.«2618E*03 

0.151676*04 

225. 

0.0359 

0.05*1 


0.07147 

0.07160 

0.77 

0.351026*03 

0.390496*03 

250. 

0.0355 

0.0356 


0.03554 

0.03556 

0.10 

0.3°465E*02 

0.354656*02 

275. 

0.035C 

0.0349 


0.00000 

0.00000 

0.00 

0.0OOO0E*00 

0.000006*00 

0 





TOTAL 

TRANSPORT t 0.2 

36«SE*0S 



0 VELOCITY CALCULATED ASSUMING A LEVEL OF NO HOTION OF Z 75. DP 
DISTANCE BETWEEN STATIONS : 3.237 KM 


88 






0 

• — 

■« GEOSTRPH1C VELOCITY 

BETWEEN STATIONS 14 AND 13 




0 2(DB) 

AV DYM( 

14) AV DYN( 13) 

SUM DYM( 14) 

SUM DYN( 13) 

VELOCITs(CM/3) 

TRANSPORT(MJ/S) 

SUM OF T ft AMI 

0 . 

0.0531 

0.0538 

0.42*22 

0-44082 

6.17 

0.35307E*04 

0.4355*6*03 

25. 

0.0^57 

0.0450 

0.38*13 

0.36*°? 

3,35 

0.237156*04 

-0.30*526*04 

50. 

0.0418 

0.0423 

0.34041 

0.34120 

3.04 

0.158°0E*Q4 

-0.54667E*04 

75. 

0 . 0:07 

0.0403 

0.2*862 

0.2°8°4 

1.24 

0 . 1 S 02 °E *02 

-0.705566*04 

100 . 

0.0186 

0 0387 

0.258°! 

0.2585° 

-1 .20 

-0.113156*04 

-0.707076*04 

125. 

0.0378 

0.0378 

0 . 22*11 

0.21 ° 8 S 

-1.85 

-0.1467*6*04 

-0.593?2E*04 

150. 

0.0373 

0.0373 

0.18258 

0.18203 

- 2.11 

-0.153816*04 

-0.447136*04 

175. 

0.0370 

0.0368 

0.14523 

0.14470 

-2.04 

-0.1J43°E* 04 

-0.2?330E*04 

200 . 

0.03*6 

0.03*4 

0.10827 

0.10786 

-1.58 

-0.95443E* 03 

-0.158?1E*04 

r:s. 

0.03b 1 

0.03*0 

0.0716° 

0.07144 

-0 °° 

-0.500 78E*03 

-0.634666*03 

250. 

0.0356 

0.0355 

0.03556 

0.03547 

-0.36 

-0.1338SE*03 

-0.133886*03 

275. 

0.034* 

0.0348 

0.00000 

0.00000 

0 00 

0.00000E* 00 

0.OOOOOE*OQ 

0 




TOTAL 

TRANSPORT « 0. 4355°E*0I 



0 VELOCITY CALCULATED ASSUMING a LEVEL OF no MOTION of 275. D 8 
DISTANCE BETWEEN STATIONS : 2.°67 KM 


0 CEOSTRPHIC VELOCITY BETWEEN STATIONS 13 AND 12 m.m 

0 2(DP> AV DVM( 13) AV DVN( 12) SUM D/N( 13) SUM DYNC 12) VELOCITY(CH/S) TRANSPORT(M3/3) SUM OF TRANS. 


0 . 

C. 05 58 

0.C55* 

0.U18? 

0 .]448 Q 

7.36 

0.60 *8 IE *04 

0.°0320E* 04 

25. 

0 0453 

0.C466 

0.08805 

0.08*28 

3.00 

0.235* 7E *04 

0,2*5 3°E* 04 

55. 

C.C423 

0 . 0 * 2 ? 

0.04225 

0 04267 

1.02 

0 S°72CE*03 

0.5°720E*03 

75. 

0.0403 

0.0301 

0.00000 

0.00000 

0 00 

0.CCOOOE * 00 ' 

0. OOOOCE* 00 





TOTAL 

TRANSPORT s 

0.903206*04 



0 VELOCITY CALCULATED ASSUMING A LEVEL OF NO MOTION OF 75. DB 
DISTANCE BETWEEN STATION’S * 4 ,6 °3 KM 


0 GECSTRFMIC VELOCITY BETWEEN STATIONS 12 AND 11 


0 203) 

AV O/NC 12) 

AV DVM( 11) SUM 0VM( 12) 

S’JM 0YN< 11) 

VELOCITV(CM/S 

) TRANSPORT(MI/S) 

SUM OF TRANS 

0 . 

0.C5S* 

0.05*2 0.1448° 

0-14612 

3.07 

C.26°02E*04 

0.365S2E*04 

25 

0.04** 

C. 0 *« "2 0.08°28 

0.08°°' 

1.62 

0.°4522E*03 

0.®5°°5E*03 

50. 

0.0427 

0.C427 0.042*7 

0.04268 

0.03 

0.147386*02 

0.14738E*02 

75. 

0.C3C1 

0.0303 0.00000 

0.00000 

0.00 

O.OOOOOE*OD 

0.000006*00 

C 



total 

TRANSPORT * 0. 

36582E*04 


0 VSL 

OCiTY CALCULA 

TED ASSUMING A LEVEL OF NO MOTION OF 75. DB 




c;: 

;TANCE BETWEEN 

STATIONS : 4.603 KM 






C 


•”»» GEC 

istdpmic velocity between s 

TAT 

IONS 11 AND 

10 




0 21 

IDP) 

AV DYNC ID 

AV DYNC 10) SUM DVN( 

ID 

SUM DYN( 

1 C) 

VELOCITV(CM/S) 

TRANSPORT(K3/S) 

SUM OF TRANS 


0. 

0.05*2 

0.0535 0. 14612' 


0.13°*° 


-15.8? 

-O.I45°OE*CS 

-0.2342?E*05 


25. 

0 04 72 

0.0447 0. C8°®3 


0.0861° 


-° .23 

-0.70°°7E*04 

-0.883866*04 


e C. 

0.0427 

0.0415 0.04268 


0.04147 


- 2 .°° 

-0.1738°E*C4 

-0.1738°E*04 


75. 

C.0303 

0.0141 O.OOOVO 


0.00000 


0 00 

0.00000E *00 

0.00000E*00 

Q 





TOTAL 

TRANSPORT * -0.2 

342*6*05 


c 

velocity CALCULAT 

ED ASSUMING A LEVEL t r NO 

MOT 

ION OF 75. 

DB 





DI 

STANCE BE TWEE’) 

STATIONS • 4.646 KM 













APPENDIX B. MEAN ACOUSTIC DOPPLER CURRENT PROFILER 

(ADCP) VELOCITIES 


ADCP STATISTICS 


(Note: Negative velocities are into the Bay, 

positive velocities are out of the Bay. ) 

GROUP: 1 

MEAN DISTANCE FROM NORTH (KM): A. 313 



MEAN 

STANDARD 

STANDARD 

DF,PTH(M) 

VELOCITY(CM/S) 

DEVIATI0N(CM/S) 

ERROR(CM/S) 

22 

-2.26 

10.54 

4. 71 

38 

-0.95 7.94 

GROUP: 2 

MEAN DISTANCE FROM NORTH (KM): 8.620 

3.55 


MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION CM/S) 

ERROR(CM/S) 

22 

-4.48 

10. 99 

4.92 

38 

-3.02 

10. 84 

4. 85 

54 

-1.29 7.84 

GROUP: 3 

MEAN DISTANCE FROM NORTH (KM): 12.950 

3.50 


MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION CM/S) 

ERROR(CM/S) 

22 

-9.46 

10. 67 

4. 77 

38 

-6.40 

9. 40 

4.20 

54 

-5.94 

9. 57 

4. 28 

70 

-4.71 5.60 

GROUP: 4 

MEAN DISTANCE FROM NORTH (KM): 17.324 

2.50 


MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION CM/S) 

ERROR(CM/S) 

22 

-14.97 

3. 86 

1. 72 

38 

-13.22 

3. 22 

1.44 

54 

-10.58 

5. 02 

2. 24 


90 





70 

-8.57 

5.41 

2.42 

86 

-4. 28 

4. 97 

2. 22 

102 

-2. 44 

3. 94 

1. 76 

118 

-0. 86 

5.45 

2.44 

134 

1. 03 

6. 66 

2. : 98 

150 

1. 64 

7. 57 

3. 39 

166 

1. 94 

6. 89 

3. 08 

182 

2. 10 

7. 11 

3. 18 

198 

1. 78 

6. 70 

3. 00 

214 

1.31 

6. 07 

2. 71 

230 

0. 63 

5. 23 

2. 34 

246 

-0. 81 

4. 70 

2. 10 

262 

0.07 

4. 78 

2. 14 

278 

0. 70 

5. 29 

2. 36 

294 

-0. 10 

4. 67 

2. 09 

310 

-1. 92 

3. 50 

1. 56 

326 

-2. 97 

2. 82 

1,26 

342 

-2.96 

2. 01 

0. 90 

358 

-3. 24 

1. 45 

0. 65 

374 

-4. 27 

1. 40 

0. 63 

390 

-4.52 

1. 99 

0. 89 

406 

-4. 62 

2. 36 

1. 06 

422 

-4. 64 

3.42 

1.53 



GROUP: 5 




MEAN DISTANCE FROM 

NORTH (KM): 21.809 



MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITYCCM/S) 

DEVIATION CM/S) 

ERROR(CM/S) 

22 

-11.47 

3. 81 

1. 70 

38 

-9. 81 

3. 09 

1. 38 

54 

-8. 36 

4. 50 

2. 01 

70 

-6.42 

2. 36 

1. 05 

86 

-4. 99 

2,87 

1. 28 

102 

-8. 37 

3. 62 

1. 62 

118 

-8. 86 

5. 29 

2. 36 

134 

-8. 47 

6. 76 

3. 02 

150 

-6. 66 

5. 64 

2.52 

166 

-4,24 

3. 23 

1.45 

182 

-3. 19 

3.59 

1. 60 

198 

-1. 99 

3. 66 

1.64 

214 

-1. 17 

2. 85 

1. 27 

230 

-1.48 

2. 54 

1. 14 

246 

-1. 71 

3,15 

1.41 

262 

-2,06 

3. 79 

1. 69 

278 

-0. 63 

3. 45 

1.54 

294 

-0,59 

3. 69 

1. 65 

310 

-1.59 

3. 19 

1.42 

326 

-2.03 

2.42 

1.08 

342 

-2. 11 

2. 16 

0. 97 

358 

-2. 60 

2. 44 

1. 09 

374 

-2. 30 

1. 93 

0. 86 

390 

-2. 03 

2. 56 

1. 14 

406 

-2. 24 

4. 21 

1. 88 


91 


422 


-3.65 


3. 89 


1. 74 


GROUP: 6 

MEAN DISTANCE FROM NORTH (KM): 25.809 



MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION CM/S) 

ERROR(CM/S) 

22 

-8. 78 

5. 77 

2.58 

38 

-8. 19 

4. 36 

1. 95 

54 

-6.91 

4. 76 

2. 13 

70 

-6.35 

4. 97 

2. 22 

86 

-6. 62 

5.45 

2.44 

102 

-6.49 

4. 20 

1.88 

118 

-3.62 

2. 67 

1. 19 

134 

-2. 34 

3.54 

1.59 

150 

-1.25 

4.87 

2. 18 

166 

-1.31 

6.53 

2. 92 

182 

-0.84 

7. 10 

3. 17 

198 

-0.48 

6. 26 

2. 80 

214 

-0. 74 

6. 84 

3. 06 

230 

0. 10 

7. 37 

3. 30 

246 

-0. 70 

6. 68 

2 99 

262 

-3. 02 

7.26 

3. 25 

278 

-3. 34 

6. 10 

2. 73 

294 

-2. 71 

4. 39 

1.96 

310 

-1.91 

3. 84 

1. 72 

326 

-2.56 

3.31 

1.48 

342 

-2. 69 

2.47 

1. 10 

358 

-2. 76 

3. 62 

1. 62 

374 

-2.56 

4. 07 

1. 82 

390 

-2. 61 

4. 28 

1.91 

406 

-2. 10 

4. 06 

1. 82 

422 

-1. 14 

3. 19 

1.42 



GROUP: 7 

MEAN DISTANCE FROM 

NORTH (KM): 29.986 



MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION(CM/S) 

ERROR(CM/S) 

22 

5.46 

3.46 

1.55 

38 

5. 07 

6. 89 

3.08 

54 

1. 79 

5.97 

2. 67 

70 

-1.94 

4. 89 

2. 19 

86 

-5. 37 

5. 89 

2. 63 

102 

-10.12 

7. 91 

3.54 

118 

-13.71 

9. 17 

4. 10 

134 

-12.88 

6. 73 

3. 01 

150 

-11. 07 

6. 02 

2.69 

166 

-10.56 

5. 84 

2. 61 

182 

-10.25 

4.53 

2. 03 

198 

-8. 94 

4. 05 

1. 81 

214 

-8.69 

3.50 

1. 57 


92 





230 

-8. 96 

4. 15 

1. 86 

246 

-8. 77 

4.43 

1. 98 

262 

-8. 17 

3. 29 

1.47 

278 

-8. 14 

2. 88 

1. 29 

294 

-8. 59 

2. 34 

1. 05 

310 

-7. 85 

1. 62 

0. 72 

326 

-6. 83 

2. 34 

1.04 

342 

-6. 00 

3. 68 

1.64 

358 

-5. 20 

4. 11 

1. 84 

374 

-5.51 

3. 39 

1.52 

390 

-6. 21 

2. 70 

1. 21 

406 

-5. 83 

2. 51 

1. 12 

422 

-5.21 2.87 

GROUP: 8 

MEAN DISTANCE FROM NORTH (KM): 34.476 

1.28 


MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATIONCM/S) 

ERROR(CM/S) 

22 

10.48 

7. 81 

3.49 

38 

6. 17 

6. 14 

2. 75 

54 

3. 32 

6. 80 

3.04 

70 

3. 92 

7. 81 

3.49 

86 

-1.56 

8. 79 

3.93 

102 

-4. 17 

10.50 

4. 70 

118 

-4. 34 

10. 49 

4.69 

134 

-2.65 

9.28 

4. 15 

150 

-1.89 

9. 24 

4. 13 

166 

-2.09 

9.52 

4.26 

182 

-2. 35 

8. 21 

3. 67 

198 

-2. 36 

6.45 

2.89 

214 

-1. 00 

4. 76 

2. 13 

230 

0. 03 

3.53 

1.58 

246 

-0. 98 

3. 66 

1. 64 

262 

-1. 08 

3. 77 

1. 69 

278 

-1.41 

3. 41 

1.53 

294 

-1. 14 

2. 61 

1. 17 

310 

-1. 28 

1. 97 

0. 88 

326 

-1. 10 

1. 37 

0. 61 

342 

-1.44 

2. 25 

1. 00 

358 

-1. 78 

2. 84 

1. 27 

374 

-2. 08 

2. 83 

1. 27 

390 

-1. 89 

3. 13 

1.40 

406 

-1. 04 

3. 49 

1. 56 

422 

-0. 56 

GROUP: 9 

MEAN DISTANCE FROM 

3. 36 

NORTH (KM): 38.597 

1. 50 


MEAN 

STANDARD 

STANDARD 

DEPTH(H) 

VELOCITY(CM/S) 

DEVIATION CM/S) 

ERROR(CM/S) 

22 

7. 02 

7. 88 

3.52 


93 


38 

-0. 32 

6. 39 

2. 86 

54 

-0. 98 

5. 19 

2. 32 

70 

-3. 28 

7.51 

3. 36 


94 


APPENDIX C. MONTEREY BAY MODEL BOUNDARY CONDITIONS 


MEAN (UNCORRECTED) VELOCITIES ALONG MOUTH SECTION 


POINT 

DIST(KM) FROM SOUTH 

VI 

V2 

DEPTH(M) 


0.33573E+02 

0.16357E+01 

0.OOOOOE+OO 

72 

2 

0.29217E+02 

0.14175E+01 

0.00000E+00 

91 

3 

0.24861E+02 

0.12107E+02 

0.00000E+00 

165 

4 

0.21431E+02 

0.12209E+02 

0,39263E+01 

483 

5 

0.19022E+02 

0.12831E+01 

-0.46778E+01 

971 

6 

0.14852E+02 

0.84049E+01 

0. 78001E+00 

1528 

7 

0.99097E+01 

-0.54945E+01 

0. 34386E+01 

1501 

8 

0.65941E+01 

-0.64503E+01 

0. 12990E+01 

940 

9 

0.26277E+01 

-0.11061E+02 

0. OOOOOE+OO 

413 



MEAN V VELOCITIES PROJECTED (INTERPOLATED) ON MODEL GRID 

GRID 

POINT 

DIST(KM) FROM SOUTH 

VI 

V2 

MODEL DEPTH(M) 

2 

0.28500E+02 

0.14175E+01 

0. OOOOOE+OO 

90 

3 

0. 21000E+02 

0.12107E+02 

0. OOOOOE+OO 

162 

4 

0. 19000E+02 

0.12209E+02 

0. 39263E+01 

450 

5 

0.12500E+02 

0.12831E+01 

-0. 46778E+01 

900 

6 

0.11500E+02 

0.84049E+01 

0., 78001E+00 

1080 

7 

0. 75000E+01 

-0.54945E+01 

0. 34386E+01 

900 

8 

0.65000E+01 

-0.64503E+01 

0. 12990E+01 

540 

9 

0.20000E+01 

-0.11Q61E+02 

0. OOOOOE+OO 

157 


POINT 

MODEL BOUNDARY U (INTERPOLATED) 
DIST(KM) FROM SOUTH VI 

VELOCITIES 

V2 

65 

0.32000E+02 

0.. 59060E+00 

0.OOOOOE+OO 

64 

0.31500E+02 

0.70873E+00 

0.OOOOOE+OO 

63 

0.31000E+02 

0.82685E+00 

0.OOOOOE+OO 

62 

0.30500E+02 

0.94497E+00 

0.OOOOOE+OO 

61 

0.30000E+02 

0.10631E+01 

0.OOOOOE+OO 

60 

0.29500E+02 

0.11812E+01 

0.OOOOOE+OO 

59 

0.29000E+02 

0. 12993E+01 

0.00000E400 

58 

0.28500E+02 

0.14175E+01 

0.OOOOOE+OO 

57 

0,28000E+02 

0.21301E+01 

0.OOOOOE+OO 

56 

0.27500E+02 

0.28427E+01 

0.OOOOOE+OO 

55 

0.27000E+02 

0.. 35554E+01 

0.OOOOOE+OO 

54 

0.26500E+02 

0.42680E+01 

0.OOOOOE+OO 

53 

0., 26000E+02 

0.49806E+01 

0.OOOOOE+OO 

52 

0.25500E+02 

0.56933E+01 

0.OOOOOE+OO 

51 

0.25000E+02 

0.64059E+01 

0.OOOOOE+OO 

50 

0.. 24500E+02 

0.71185E+01 

0.OOOOOE+OO 

49 

0.; 24000E+02 

0. 78312E+01 

0.OOOOOE+OO 

48 

0., 23500E+02 

0.85438E+01 

0.OOOOOE+OO 

47 

0.23000E+02 

0.92565E+01 

0.OOOOOE+OO 

46 

0.. 22500E+02 

0.99691E+01 

0.OOOOOE+OO 
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45 

0.22000E+02 

0. 10682E+02 

44 

0.21500E+02 

0. 11394E+02 

43 

0.21000E+02 

0. 12107E+02 

42 

0.20500E+02 

0. 12132E+02 

41 

0.20000E+02 

0. 12158E+02 

40 

0.19500E+02 

0.12183E+02 

39 

0.19000E+02 

0.12209E+02 

38 

0.18500E+02 

0. 11368E+02 

37 

0.18000E+02 

0.10528E+02 

36 

0.17500E+02 

0.96874E+01 

35 

0.17000E+02 

0.88470E+01 

34 

0.16500E+02 

0.80065E+01 

33 

0.16000E+02 

0.71661E+01 

32 

0.15500E+02 

0.63257E+01 

31 

0.15000E+02 

0.54853E+01 

30 

0.14500E+02 

0.46448E+01 

29 

0.14000E+02 

0.38044E+01 

28 

0.13500E+02 

0.29640E+01 

27 

0.13000E+02 

0.21236E+01 

26 

0.12500E+02 

0.12831E+01 

25 

0.12000E+02 

0.48440E+01 

24 

0.11500E+02 

0.84049E+01 

23 

0.11000E+02 

0.66675E+01 

22 

0.10500E+02 

0.49301E+01 

21 

0.10000E+02 

0.31926E+01 

20 

0.95000E+01 

0.14552E+01 

19 

0.90000E+01 

-0.28223E+00 

18 

0.85000E+01 

-0.20197E+01 

17 

0.80000E+01 

-0.37571E+01 

16 

0.75000E+0 1 

-0.54945E+01 

15 

0.70000E+01 

-0.59724E+01 

14 

0.65000E+01 

-0.64503E+01 

13 

0.60000E+01 

-0.69626E+01 

12 

0. 55000E+01 

-0.74750E+01 

11 

0.50000E+01 

-0.79873E+01 

10 

0.45000E+01 

-0.84996E+01 

9 

0.40000E+01 

-0.90119E+01 

8 

0. 35000E+01 

-0.95243E+01 

7 

0.30000E+01 

-0.10037E+02 

6 

0.25000E+01 

-0.10549E+02 

5 

0.20000E+01 

-0.11061E+Q2 

4 

0.15000E+01 

-0.82959E+01 

3 

0.10000E+01 

-0.55306E+01 

2 

0. 50000E+00 

-0.27653E+01 

1 

0.OOOOOE+OO 

0.OOOOOE+OO 


0. 00000E+00 
0. 00000E+00 
0. OOOOOE+OO 
0. 98157E+00 
0. 19631E+01 
0. 29447E+01 
0. 39263E+01 
0. 32644E+01 
0. 26026E+01 
0. 19407E+01 
0. 12789E+01 
0. 61703E+00 
-0. 44824E-01 
-0. 70668E+00 
-0. 13685E+01 
-0.20304E+01 
-0.26922E+01 
-0.33541E+01 
-0.40159E+01 
-0.46778E+01 
-0.19489E+01 
0.78001E+00 
0.11123E+01 
0.14446E+01 
0.17770E+01 
0.21093E+01 
0.24416E+01 
0.27739E+01 
0.31063E+01 
0.34386E+01 
0.23688E+01 
0.12990E+01 
0.11547E+01 
0.10103E+01 
0.86600E+00 
0.72167E+00 
0.57733E+00 
0.43300E+00 
0.28867E+00 
0.14433E+00 
0.00000E+00 
0.OOOOOE+OO 
0.OOOOOE+OO 
0.OOOOOE+OO 
0.OOOOOE+OO 
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DERIVATION OF MASS BALANCED FLOW FIELD 


RAW VALUES MASS BALANCED VALUES 

(M**2/S) (M**2/S) 


GRID POINT 

VI 

V2 

VI 

V2 

1 

0.0 

0.0 

0. 0 

0.0 

3 

0.5 

0. 0 

0. 3 

0.0 

5 

0. 6 

0. 0 

0.4 

0. 0 

7 

0.7 

0.0 

0. 5 

0. 0 

9 

1.4 

0.0 

0. 9 

0. 0 

11 

2. 1 

0. 0 

1.4 

0. 0 

13 

2. 8 

0. 0 

1. 8 

0. 0 

15 

3.6 

0. 0 

2.4 

0. 0 

17 

4.3 

0. 0 

2. 8 

0. 0 

19 

5.0 

0. 0 

3. 3 

0. 0 

21 

5. 7 

0. c 

3. 7 

0. 0 

23 

6. 1 

0. 0 

4. 0 

0. 0 

25 

6. 1 

4. 9 

4. 0 

3. 6 

27 

5. 7 

11. 0 

3. 7 

8. 1 

29 

4. 8 

9.6 

3. 1 

7. 1 

31 

4. 0 

3. 9 

2. 6 

2.9 

33 

3. 2 

-5.4 

2. 1 

-7.9 

35 

2. 3 

-16. 0 

1.5 

-23.4 

37 

1.5 

-26. 0 

1. 0 

-38.0 

39 

0.6 

-40. 0 

0.4 

-58.4 

41 

4.2 

7.4 

2.8 

5.5 

43 

2.5 

16.0 

1.6 

11. 8 

45 

0. 7 

26. 0 

0. 5 

19. 2 

47 

-1.0 

34. 0 

-1.5 

25. 2 

49 

-2. 7 

40. 0 

-4. 2 

29. 6 

51 

-3. 2 

11. 0 

-4. 9 

8. 1 

53 

-3. 7 

5. 1 

-5. 7 

3. 8 

55 

-4. 2 

2. 5 

-6.5 

1. 8 

57 

-4. 8 

0. 9 

-7.4 

0. 6 

59 

-5. 3 

0. 3 

-8. 2 

0. 2 

61 

-4. 1 

0. 0 

-6. 3 

0. 0 

63 

0. 0 

0. 0 

0. 0 

0. 0 

NET 

TRANSPORT 

39.4 

85. 1 

0. 0 

0. 0 


LINEAR CORRECTION COEFFICIENTS APPLIED TO: 

NEGATIVE VI VALUES = 1.5450 
POSITIVE VI VALUES = 0.6550 
NEGATIVE V2 VALUES = 1.4605 
POSITIVE V2 VALUES = 0.7399 


97 


MODEL BATHYMETRY ALONG OPEN BOUNDARY 


DISTANCE(KM) FROM SOUTH MODEL DEPTH (M) 


31.5 

90. 0 

31.0 

90. 0 

30.5 

90. 0 

30. 0 

90. 0 

29.5 

90. 0 

29. 0 

90. 0 

28.5 

90. 0 

28. 0 

90. 0 

27.5 

90. 0 

27. 0 

90. 0 

26.5 

90. 0 

26. 0 

90.0 

25.5 

90.0 

25.0 

90. 0 

24.5 

90.0 

24.0 

94.5 

23.5 

99. 0 

23. 0 

103.5 

22.5 

108.0 

22. 0 

126.0 

21. 5 

144. 0 

21. 0 

162. 0 

20.5 

180. 0 

20. 0 

270. 0 

19. 5 

360. 0 

19. 0 

450. 0 

18.5 

540. 0 

18. 0 

585.0 

17.5 

630. 0 

17. 0 

675.0 

16.5 

720. 0 

16. 0 

720.0 

15.5 

720.0 

15.0 

720. 0 

14.5 

720.0 

14. 0 

765.0 

13.5 

810. 0 

13. 0 

855. 0 

12.5 

900. 0 

12. 0 

990. 0 

11. 5 

1080.0 

11. 0 

1170.0 

10.5 

1260.0 

10. 0 

1260.0 

9.5 

1260.0 

9. 0 

1260.0 

8.5 

1260.0 

8. 0 

1080.0 

7.5 

900. 0 

7.0 

720.0 

6.5 

540. 0 

6.0 

461.2 
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5.5 

382.5 

5.0 

303. 7 

4.5 

225.0 

4. 0 

213. 7 

3.5 

202.5 

3. 0 

191. 2 

2.5 

180. 0 

2. 0 

157.5 

1.5 

135. 0 

1 . 0 

112.5 

0.5 

90. 0 
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APPENDIX D. MEAN BRUNT-VAISALA FREQUENCIES AND ROSSBY 

RADII 


STATION: 1 LAT: 36 55.2 LON: 122 5.8 

BOTTOM DEPTH = 45.5 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.1174E-01 
EXTERNAL ROSSBY RADIUS (KM) = 241.85 

FIRST INTERNAL ROSSBY RADIUS (KM) = 1. 95 


STATION: 2 LAT: 36 52.8 LON: 122 5.9 

BOTTOM DEPTH = 75. 2 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.1063E-01 
EXTERNAL ROSSBY RADIUS (KM) = 311.12 

FIRST INTERNAL ROSSBY RADIUS (KM) = 2.92 


STATION: 3 LAT: 36 50.5 LON: 122 5.9 

BOTTOM DEPTH = 91. 1 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.9942E-02 
EXTERNAL ROSSBY RADIUS (KM) = 342.38 

FIRST INTERNAL ROSSBY RADIUS (KM) = 3.30 


STATION: 4 LAT: 36 48. 1 LON: 122 5.8 

BOTTOM DEPTH = 186. 3 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.7315E-02 
EXTERNAL ROSSBY RADIUS (KM) = 489.65 
FIRST INTERNAL ROSSBY RADIUS (KM) = 4.97 


STATION: 5 LAT: 36 46.8 LON: 122 5.8 

BOTTOM DEPTH = 669. 8 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.4895E-02 
EXTERNAL ROSSBY RADIUS (KM) = 928.38 

FIRST INTERNAL ROSSBY RADIUS (KM) = 11.96 


STATION: 6 LAT: 36 45.5 LON: 122 5.8 

BOTTOM DEPTH = 1023.6 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.4220E-02 
EXTERNAL ROSSBY RADIUS (KM) = 1147.67 
FIRST INTERNAL ROSSBY RADIUS (KM) = 15.76 
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STATION: 7 LAT: 36 42.3 LON: 122 5.8 

BOTTOM DEPTH = 1582.0 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.3597E-02 
EXTERNAL ROSSBY RADIUS (KM) = 1426.77 
FIRST INTERNAL ROSSBY RADIUS (KM) = 20.76 


STATION: 8 LAT: 36 41.1 LON: 122 3.6 

BOTTOM DEPTH = 891.3 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.4457E-02 
EXTERNAL ROSSBY RADIUS (KM) = 1070.94 
FIRST INTERNAL ROSSBY RADIUS (KM) = 14.49 


STATION: 9 LAT: 36 40.3 LON: 122 2. 1 

BOTTOM DEPTH = 594. 6 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.5098E-02 
EXTERNAL ROSSBY RADIUS (KM) = 874.71 

FIRST INTERNAL ROSSBY RADIUS (KM) = 11.06 


STATION: 10 LAT: 36 38.8 LON: 121 59. 1 
BOTTOM DEPTH = 83. 2 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.1041E-01 
EXTERNAL ROSSBY RADIUS (KM) = 327.13 

FIRST INTERNAL ROSSBY RADIUS (KM) = 3.16 


STATION: 11 LAT: 36 41.2 LON: 121 58.2 
BOTTOM DEPTH = 93. 1 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.1005E-01 
EXTERNAL ROSSBY RADIUS (KM) = 346.09 
FIRST INTERNAL ROSSBY RADiUS (KM) = 3.41 


STATION: 12 LAT: 36 43.6 LON: 121 57.4 
BOTKM DEPTH = 93. 1 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.9903E-02 
EXTERNAL ROSSBY RADIUS (KM) = 346.09 

FIRST INTERNAL ROSSBY RADIUS (KM) = 3.36 


STATION: 13 LAT: 36 46. 1 LON: 121 57. 9 
BOTTOM DEPTH = 648. 2 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.4821E-02 
EXTERNAL ROSSBY RADIUS (KM) = 913.24 

FIRST INTERNAL ROSSBY RADIUS (KM) = 11.40 
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STATION: 14 LAT: 36 47.6 LON: 121 58.6 
BOTTOM DEPTH = 297. 5 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.5949E-02 
EXTERNAL ROSSBY RADIUS (KM) = 618.68 

FIRST INTERNAL ROSSBY RADIUS (KM) = 6.46 


STATION: 15 LAT: 36 49.3 LON: 121 59. 1 
BOTTOM DEPTH = 428. 2 

MEAN BRUNT-VAISALA FREQ (CPS'* = 0. 5447E-02 
EXTERNAL ROSSBY RADIUS (KM) = 742.25 

FIRST INTERNAL ROSSBY RADIUS (KM) = 8.51 


STATION: 16 LAT: 36 51.9 LON: 121 59.9 
BOTTOM DEPTH = 59.4 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.1177E-01 
EXTERNAL ROSSBY RADIUS (KM) = 276.58 

FIRST INTERNAL ROSSBY RADIUS (KM) = 2. 55 


STATION: 17 LAT: 36 54.1 LON: 122 0.7 

BOTTOM DEPTH = 27. 8 

MEAN BRUNT-VAISALA FREQ (CPS) = 0.1490E-01 
EXTERNAL ROSSBY RADIUS (KM) = 189.10 

FIRST INTERNAL ROSSBY RADIUS (KM) = 1.51 
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